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The coherent and collective nature of a Bose-Einstein condensate can enhance or suppress physical 
processes. Bosonic stimulation enhances scattering in already occupied states which leads to matter 
wave amplification, and the suppression of dissipation leads to superfluidity. In this summer school 
notes we present several experiments where enhancement and suppression have been observed and 
discuss the common roots of and differences between these phenomena. 
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When a gas of bosonic atoms is cooled below the transition temperature of Bose-Einstein condensation, it profoundly 
changes its properties. The appearance of a macroscopically occupied quantum state leads to a variety of new 
phenomena which set quantum fluids apart from all other substances. Fritz London even called them the fourth state 
of matter p[ . 

Many of the key concepts in quantum fluids were derived from studying the weakly interacting Bose gas, for which 
rigorous theoretical treatments were possible @,g. In 1995, with the discovery of BEC in a dilute gas of alkali 
atoms it became possible to study such a system experimentally. Since then, the interplay between theory and 

experiment has been remarkable jS. 

This paper is the fourth major review paper of our group which describes the new techniques and the physics of 
Bose-Einstein condensation p|-|io|| . These review papers together give a comprehensive overview of the topics to which 
our group has contributed. Each contribution is self-contained, but we avoided major overlap with previous review 
papers. The topic of these Cargese summer school notes is enhancement and suppression of physical processes in a 
Bose-Einstein condensate. 

Many phenomena in Bose condensates involve such enhancement or suppression. Our recent experiments include 
the enhancement and suppression of elastic collisions of impurity atoms 11 1, the suppression of dissipation due to 
superfluidity JlJ,[l3| , and the enhancement [Q and suppression jl5| of light scattering. The common discussion of 
these phenomena leads to a better understanding of the underlying principles. We draw analogies between light 
scattering and particle scattering, between microscopic and macroscopic superfluidity. We show that a condensate 
responds very differently to two different ways of momentum transfer, light scattering and spontaneous emission [fl6|| . 
We discuss light scattering in both the linear and nonlinear regime where bosonically enhanced Rayleigh scattering 
led to the amplification of either atoms jlTj or light (l8) and work out the relationship between these two processes. 
Finally, the section on matter wave amplification of fermions discusses the relevance of symmetry and long coherence 
time and its relation to quantum statistics. 



I. SCATTERING OF LIGHT AND MASSIVE PARTICLES 



Before we discuss light scattering and collisions in a BEC, we want to derive some simple general expressions based 
on Fermi's golden rule. This will help to see the similarities and differences between the various processes. When a 
condensate scatters a photon or material particle, the scattering is described by the Hamiltonian 

H' = C H^nCka m 8i +n -k-m- (1) 

k,l,m,n 

Here c k (c k ) is the destruction (creation) operator for the scattered particles (which can be photons or massive 
particles), and a k (at) is the destruction (creation) operator for atomic plane waves of wavevector k (free particle 
states). The strength of the coupling is parametrized by the coefficient C (which in general may depend on the 
momentum transfer), and the 5 function guarantees momentum conservation (with an uncertainty of h/D, where D 
is the dimension of the condensate). 

We consider the scattering process where a system with Nq atoms in the condensate ground state and N q quasi- 
particles with wavevector q scatters particles with incident wavevector k into a state with wavevector k — q. The 
initial and final states are[] 

\i) = \n k ,n k - q ;N a ,N q ) 

\.f) = K - 1, ri k - q + 1; JV - 1, N q + 1) (2) 

respectively, where n k denotes the population of scattering particles with wavevector k. It should be emphasized that, 
due to the interatomic interactions, the quasi-particles with occupation N q are not the plane waves created by the 
operator a' q , but the quanta of collective excitations with wavevector q. 
The square of the matrix element M± between the initial and final state is 



1 This choice of final states implies that we neglect scattering between quasi-particles and consider only processes involving 
the macroscopically occupied zero-momentum state of the condensate. Formally, we replace the Hamiltonian (Eq. [j]) by 
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W+ = — \C\ 2 S(q)N [n k (N q + n k _ q + 1) - N q n k . q ]5{E k - E k _ q - huf). (7) 



|Mi| a = |</|W'|i)| a 

= |C| 2 |(7V = N — l,N q = l|p+(q)|iVo =N,N q = 0)\ 2 (N q + l)(n fc _, + l)n k (3) 

where fft(q) = ^ m o^ m + q a m is the Fourier transform of the atomic density operator at wavevector q. The static 
structure factor of the condensate is 

S(q) = ( 5 |p(q)pt(q)| 5 )/iV (4) 

where \g) = \Nq = N, N q = 0) is the BEC ground state and N is the total number of atoms. We then obtain for the 
scattering matrix element Mi 

\Mi\ 2 = \C\ 2 S(q)(N q + l)(n fc _, + l)JV n fc . (5) 

The scattering rate W\ for the process \n k , n k -q', No, N q ) — ► \n k — 1, n k - q +1;Nq — 1, N q + 1) follows from Fermi's 
golden rule as 

W x = ^\M^6{E k - E k _ q - Tiuf ) (6) 

where E k is the energy of the incident particle with wavevector k, and Uui q is the energy of quasi-particles with 
wavevector q (which we will later obtain from Bogoliubov theory). To obtain the net growth rate of N q , one has to 
include the reverse process \n k , n k ^ q ; Nq, N q ) — ► \n k + 1, n k ^ q — 1; A^o + 1, N q — 1) by which atoms scatter back into 
the condensate. The square of the matrix element Mi for this process is \C\ 2 S(q)N q n k - q {No + l)(n k + 1). The net 
rate W+ of scattering atoms from the condensate into the quasi-particle mode q is the difference of the two partial 
rates W+ = W\ — W<z. Assuming A*o 3> 1 (i.e., Nq + 1 w No), we obtain for the net rate 

For large n^. (e.g. a laser beam illuminating the condensate) the term N q n k - q can be neglected, and the dominant 
bosonic stimulation term (N q + n k —q + 1) is approximately (N + 1) with N — max(7V g , n k -q). This illustrates that 
there is no bosonic stimulation of the net rate by the least populated final state. With the dynamic structure factor 
S(q, uj) = S(q)S(uj - uf) Eq. | simplifies to 

W + = T2\°\ 2n k S (q, (E k - E k _ q )/h) N (N q + n k _ q + 1). (8) 

The rate W + in Eq. ^ is the rate for the Stokes process where E k > E k ^ q . Momentum transfer hq to the condensate 
is also possible as an anti-Stokes process where a quasi-particle with momentum — ?iq is scattered into the condensate, 
and the scattered particle gains energy. The net rate W_ for this process is obtained in an analogous way as 

W- = ^\C\ 2 S{q)N [n k {N_ q - n k . q ) - {N_ q + \)n k ^ q ]5{E k - E k _ q + huf). (9) 

The net scattering rates of Eqs. [?] and || involve the product of three terms. 

• The static structure factor S(q) represents the squared matrix element for the condensate to absorb momentum 
hq. 

• The S function denotes the density of final states. 

• The bosonic stimulation term represents stimulation by the occupation in the final state either of the scattering 
particles or the condensate. 

The interplay of these three terms is responsible for the enhancement and suppression of physical processes in a 
condensate. The properties of the condensate as an intriguing many-body system are reflected in the structure factor. 
In Sect, [n], we discuss its measurement through stimulated light scattering. The density of states is responsible for 
superfluidity because it vanishes for initial velocities of the incident particles which are smaller than the Landau 
critical velocity (Sects. IV and |v|). Finally, bosonic stimulation by the occupancy N q of final states was responsible 
for superradiance, matter wave amplification and optical amplification in a condensate (Sects. VII and VIII). 

Sects. IX and [x] broaden the above discussion. In Sect. IX we show that there is a major difference how the 
condensate affects light scattering and spontaneous emission, i.e., spontaneous emission can probe properties of the 
condensate beyond the structure factor. Eq. |?] seems to imply that enhancement of a process requires macroscopic 
population of a quantum state, i.e., bosonic quantum degeneracy. The discussion on matter wave amplification for 
fermions in Sect, [x] shows that this is not the case, and that collective enhancement is possible even for fermions when 
they are prepared in a cooperative state. 
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II. DETERMINATION OF THE STATIC STRUCTURE FACTOR BY BRAGG SPECTROSCOPY 



The matrix element or the structure factor in Eqs. M and |9| can be directly determined experimentally by light 
scattering. The density of states (the S function in Eq. 0) does not restrict the scattering since the photon energy is 
much higher than the quasi-particle energy. As a result, photons can be scattered into the full solid angle and provide 
the necessary recoil energy of the atom by a small change (« 1CP 9 ) in the frequency of the scattered photon. 

When the light scattering is stimulated by a second laser beam, the momentum transfer hq is determined by the 
angle 9 between the two laser beams: hq = 2hksm(9/2). By varying the angle, one can probe both the phonon 
and free particle regime of the Bogoliubov quasiparticles (Fig. |l|). Fig. |^ illustrates our experimental method for 
probing the response of a condensate at large momentum transfer. Counter-propagating laser beams were incident 
perpendicular to the long axis of the cigar-shaped condensate, which contained several million sodium atoms in the 
F = 1, uijp = —1 state pj| . After the momentum transfer, the condensate was allowed to freely expand, allowing 
scattered atoms to separate spatially from the unscattered atoms, as shown in the figure. 




Momentum 



FIG. 1. Energy level diagram for Bragg spectroscopy. The solid line represents the energy of an atom at momentum hq, 
which is given by the Bogoliubov dispersion relation for a homogeneous weakly- interacting Bose-Einstein condensate (Eq. jHJ). 
Excitations can be created optically by stimulated light scattering using two laser beams which are both far detuned (about 2 
GHz) from the atomic resonance. Momentum and energy are provided by absorption of one photon to a virtual excited level, 
followed by stimulated emission of a second, lower energy photon. For small momenta hq, such that hq <C Mc a the dispersion 
relation is phonon-like (linear). Here M is the mass, and c 3 is the speed of sound which is related to the interaction energy 
H by fi — Mc? a . For large momenta hq ^> Mc s the dispersion relation is particle- like (quadratic) offset from the energy of a 
free-particle excitation by a mean-field shift of (X (« 5 kHz for our experiments). 



For strong l aser p ulses and short times, the atoms undergo Rabi oscillations between the initial state and the recoil 
state (see Sec. VIII ). For longer times, we can use rate equations (Eqs. and |J), and obtain the rate W of transferring 
photons from one beam to the other as 



W/N = {W+ + W-)/N = 27r(Oii/2) 2 S{q)[S(u - u?) - S(u + u?)]. 



(10) 



The two laser beams have wavevectors k and k — q and a difference frequency u>. The two-photon Rabi frequency 
is given by — \C\ 2 n k n k ^ q (see Sec. |V| 



When uj is scanned the "spectrum" of a condensate consists of two peaks at 



(Fig. |). The strength (or 

integral) of each peak corresponds to S{q). We refer to this method as Bragg spectroscopy since the basic process is 
Bragg scattering of atoms from an optical standing wave. A full account of this method and the underlying theory 
was given in our Les Houches notes [|l0) (see also ) . 
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FIG. 2. Bragg spectroscopy. Atoms were stimulated by two counter-propagating laser beams to absorb a photon from one 
beam and emit it into the other beam, resulting in momentum transfer to the atoms, as observed in ballistic expansion after 
20 ms time of flight. The number of scattered atoms showed a narrow resonance when the difference frequency v between the 
two laser beams was varied around the recoil frequency vq of the atoms. 
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FIG. 3. Probing S(q) by Bragg spectroscopy. Momentum transfer per particle, in units of Tiq, is shown vs. the frequency 
difference ui/2tt between the two Bragg beams. The two spectra were taken at different atomic densities. The open symbols 
represent the phonon excitation spectrum for a trapped condensate (at a chemical potential jj,/h — 9.2 kHz, much larger 
than the free recoil shift of « 1.4 kHz). Closed symbols show the free-particle response of a twenty-three times more dilute 
(ballistically expanded) cloud. Lines are fits to the difference of two Gaussian line shapes representing excitation in the forward 
and backward directions. See Refs. JJ^UhJ for more details. Figure is taken from Ref. |15|. 



We observed that the scattering rate was strongly suppressed when the momentum transfer fiq became smaller 
than the speed of sound c s in the condensate (times the atomic mass M), i.e., when the light scattering excited a 
phonon and not a free particle. These observations are in agreement with the Bogoliubov theory, which obtains the 
elementary excitations or Bogoliubov quasiparticles as eigenstates of the Hamiltonian 

h 2 k 2 

H = y^^ ~^~ a l°fc + G ^ a\a\a k a m 8i +n - k ^ m . (11) 

k k,l,m,n 

The constant C can be expressed by the scattering length a as 

C = 2 -^ (12) 

where V is the condensate volume. The structure factor S(q) (Eq. |j) is the norm of the state vector p(q)^ \g) / '\TR 
with the atomic density operator p(q)^ = 'il n+q a m . Only terms involving the zero-momentum state m = yield 
significant contributions. Thus S(q) is the norm of the vector 

(alan + aXa-„)\q) . 
uJ ' ° qn ^^(al + a^)\g) = \e+) + \e-) (13) 



N 
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where we have replaced dj and do by y/N following the usual Bogoliubov formalism j| . 

To calculate the norm of |e) explicitly, we transform to Bogoliubov operators b q = u q a q + v q a~ q , where u q = cosh(j) q , 
v q = sinh^jq and tanh20 9 = p,/(hu> q +p,) where huj® = Ti 2 q 2 /2M is the free-particle kinetic energy and p is condensate's 
chemical potential. This yields S(q) = (u q — v q ) 2 . The static structure factor tends to S(q) — > hq/2mc s and vanishes 
in the long wavelength limit, as required of a zero-temperature system with finite compressibility p2| . The speed 
of sound c s is related to the chemical potential p by p = Mc 2 s . The resonances in the Bragg spectrum occur at the 
energies hui q of quasi-particles created by b q 



fuj* = ^0(^0 + 2^). (14) 

In our experiments we have observed the suppression of Bragg scattering at small angles. This implies an overall 
suppression of Rayleigh scattering which will be more pronounced when the speed of sound becomes comparable or 
larger than the recoil velocity. By integrating the static structure factor S(q) over all possible scattering angles and 
accounting for the dipolar emission pattern, we find that Rayleigh scattering from a BEC is suppressed by a factor 

pscatt _ k s , M i : _±. | / | r i 



15 fcj 9 kl 3 K_ 
8 kl + 4 k 4 L + 2 k 2 L 




tanh 1 



where fik^ is the incident photon momentum and Tik s — Mc s is the momentum of an atom moving at the speed of 
sound. 

The long-wavelength suppression of S(q) reveals a remarkable many-body effect. For free particles, the matrix 
element for momentum transfer is always 1, which reflects the fact that the operator exp(iq • r) connects an initial 
state with momentum p to a state with momentum p + ?iq with unity overlap. For an interacting Bose-Einstein 
condensate, this overlap vanishes in the long- wavelength limit. As we have discussed in previous publications p^ , |lO| |, 
this suppression of momentum transfer is due to an destructive interference between the two pathways for a condensate 
to absorb momentum ?iq and create a quasi-particle: one pathway annihilates an admixture with momentum — ?iq, 
the other creates an admixture at momentum +hq. 

The leading terms in stimulated light scattering do not depend on temperature. The rates W+ and W- (Eqs. 
and ^) are independent of the thermally excited population of quasiparticles N q and N^ q in the limit of large 
rife, nk- q 3> N± q , i.e., when the scattering is stimulated by a second laser beam. For spontaneous scattering (rik- q = 0), 
one obtains for the total scattering rate W instead of Eq. [n] 

W/N = ^\ C \ 2 S(q)n k [(1 + N q )(S(Lu - w f ) + N- q 8(u + <))] . (16) 

Absorption of and bosonic stimulation by thermally excited quasi-particles become important when the temperature 
is comparable to or larger than the quasi-particle energy hui q . Eq. |l6| is proportional to the temperature dependent 
dynamic structure factor Sr(q, u) [|3j. 



III. ATOMIC COLLISIONS IN A BOSE-EINSTEIN CONDENSATE 



Before we discuss the analogies and differences between the scattering of light and massive particles, we make 
general remarks on particle scattering. We raise the question under what conditions can two matter waves penetrate 
each other without scattering. This includes the question what happens when two atom laser beams cross each other. 

As we have seen in Sect. |, the scattering rate is proportional to the structure factor 5(q), which is a measure of 
the density-density correlations at wavevector q associated with the momentum transfer. It consists of two parts, one 
reflects the average atomic density p(q) — (g\p((l)\g) and the other one the fluctuations Sp(q) = p(q) — p(q) (see also 
Eq. H). 

S(q) = ^<s|(p(q) + 5p(q))(p*(q) + 5p\q))\g) (17) 
= ^( / 5(q)p*(q) + (5|^(q)5pt(q)| 5 )) (18) 



G 



In a system without fluctuations (<5/3(q) = 0), there is only scattering when the stationary density modulation 
allows for a momentum transfer at wavevector q ( p(q) ^ 0) and the incident particles fulfil the Bragg condition. One 
well known example is X-ray scattering off a crystal lattice. Electrons in such a lattice form stationary Bloch waves 
(superpositions of plane waves and Bragg scattered plane waves) and propagate without attenuation. Scattering only 
occurs at irregularities of the lattice or thermal fluctuations. 

A Bose-Einstein condensate might appear perfectly ordered. However, as we have seen in Sect, [n], it has density 
fluctuation similar to the classical ideal gas. The structure factor only differs from an ideal gas when the momentum 
transfer is comparable or less than the speed of sound (times the mass M). In the case of electromagnetic waves, these 
density fluctuation cause Rayleigh scattering. In close analogy, if a matter wave propagates through a condensate, 
there will be elastic scattering. On length scales larger than the healing length, the condensate has reduced density 
fluctuations, and we have seen that electromagnetic radiation of sufficiently long wavelength can propagate with only 
little scattering (Eq. |l5|). In the next section, we will find the equivalent result for long wavelength matter waves, but 
also a new effect, namely the complete suppression of scattering for velocities below the Landau critical velocity. 

The scattering between two Bose-Einstein condensates can be nicely demonstrated by creating two condensates in 
a magnetic trap, one sitting at the bottom, the other one held up in the trapping potential by a blue-detuned light 
sheet which forms a repulsive potential for atoms. This initial situation was created by splitting a condensate into two 
halves using a light sheet j24| and then shifting the center of the magnetic trap by applying a magnetic field gradient. 
When the light sheet was switched off, one condensate accelerated and slammed into the other. The violent collision 
was observed in situ by phase contrast imaging and also analyzed by absorption imaging after ballistic expansion 
(Fig. |). 
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FIG. 4. Collisions between two condensates. One condensate was held up in the magnetic trapping potential by a 
blue-detuned light sheet and suddenly released. A series of phase-contrast images (upper part) and time-of-flight absorp- 
tion images display the collision. The phase contrast pictures show the spreading of the upper condensate during acceleration, 
its collision, and merging with the stationary cloud. The time-of-flight pictures represent the velocity distribution. Before the 
scattering event, there were two distinct peaks representing the accelerated and the stationary condensate. After the collision, 
only one cloud was left which was heated up by the collision. 
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FIG. 5. Observation of elastic collisions between the condensate (lower black dot) and Bragg diffracted atoms (upper spot). 
The image is taken after a time-of-flight of 30 ms, and shows the velocity distribution after the collision. The products of the 
collision are distributed over a sphere in momentum space leading to the observed s-wave halo. The height of the image is 3.2 
mm. 



The method shown in Fig. ^ suffers from a large velocity spread of the incident condensate. A narrower velocity 
distribution can be achieved by Bragg scattering. With counter-propagating laser beams, we transferred two photon 
recoil momenta to some of the atoms. After the recoiling atoms had traversed the condensate the velocity distribution 
was probed in time-of-flight imaging and showed the characteristic s-wave halo (Fig. |). 

Such collisions between indistinguishable atoms of mass M are described by the interaction Hamiltonian in Eq. 
[Tl| . This Hamiltonian is not only responsible for the interactions between two condensates, but also for the mean- 
field energy of a single condensate (Cajajetoao), pair correlations and quantum depletion (C^2 q aJa_ g ao&o) and 
the damping of collective excitations where a quasi-particle in mode k decays into two other excitations (Beliaev 
damping, 2CJ^j a^ajakaoSi+n-k-o) or promotes an existing quasi-particle in mode I to higher energy (Landau 

damping, 4CJ2k,n 6o*n°j6*A+n-i-k)- 

Several papers discussed the outcoupling of atoms in mode k from a condensate at rest (mode 0) and only considered 
the term 2Ca k al i aoak which represents the mean-field repulsion between the two matter waves (see, for example, 
^5),^6|). This term leads to a distortion of the outcoupled beam. However, this two-mode approximation neglects the 
scattering into the empty modes, described by 

H' = 2C ^ aja^akaoSi+n-k-o (19) 

La 

In the following we want to show that this interaction leads to a scattering rate for particles in mode k 

Wtot = n av (20) 

with the condensate density no = Nq/V, the cross-section a = 87ra 2 and the relative velocity v — fik/M. Eq. 
therefore describes not only classical particles which collide with a cross-section cr, but also the quantum mechanical 
interaction between two matter waves. 

We assume that the incident velocity is much larger than the speed of sound t)>c s , Therefore, there is no distinction 
between free particles and quasi-particles. In the center-of-mass frame the relevant terms in the Hamiltonian (Eq. [Tsj ) 
are 

W = 4C]T' alal q a k/2 a_ k/2l (21) 



where the prime indicates summation over one hemisphere. The additional factor of 2 (compared to Eq. |l9| ) accounts 
for the fact that the final state (q, — q) appears twice in the summation (as q, — q and — q, q), we only count the number 
of states in one hemisphere). Considering the collision of one particle with momentum Tik with a pure condensate, 
we write the initial and final states \Nk/2, N q ; N_ k /2, N_ q ) (see Eq. ||) as 

\i) = |1,0;JV ,0) 

|/) = |0,1;JV - 1,1). (22) 

The square of the matrix element for the Hamiltonian (Eq. |2l|) between those two states is \M\ 2 = 16\C\ 2 Nq. From 
Eq. H (Fermi's golden rule) the scattering rate follows as 

W tot = -y2\M\ 2 S(E f -E i ) 



f 



(23) 



with the final and initial energies £/=2x h 2 q 2 /2M, Ei — 2 x h 2 (k/2) 2 /2M. The density of final states dp/dE/ is 
equal to half the density of single particle states evaluated at -Esingic = Ef/2. The density of single particle states at 
energy E is 
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dp 1 V /->W\ 3 /2 

X 



dE 2 (2tt) 3 V 

The additional factor of (1/2) comes from the fact that for pairs of indistinguishable particles, we only count the 
number of states in one hemisphere. 
Thus, the density of final states is 



dp 1 dp 



dE f 2 dE 



E=E S ,2 (2^) 3 2 V h 2 ) 



(£//2) V2 (25) 



Substituting Eq. |2^ into Eq. ^3] gives q — k/2 and finally Eq. 

In the center-of-mass frame, the scattered particles occupy a shell in momentum space at q = k/2 (see Fig. ||). The 
quantum-mechanical origin of the relative velocity v in Eq. ^o] is the density of final states which is proportional to 
\[E oc v. The scattering of atoms into empty modes is not described by the Gross-Pitaevskii equation (which only 
describes macroscopically occupied modes) , but can be accounted for by introducing a complex scattering length into 
the Gross-Pitaevskii equation [^7| . 

For the scattering of impurity atoms with the same mass as the condensate atoms the constant C in Eq. |l| is 

C=^, (26) 
or more generally C = 2nh 2 a/ pV where p is the reduced mass and a is now the scattering length for collisions between 



condensate and impurity atoms. The factor of two difference between Eqs. |26| and 12 is necessary to avoid double 
counting of identical atom pairs. If we repeat the above derivation for impurity atoms we obtain a cross-section of 
a = Aira 2 , whereas the cross-section for indistinguishable particles is twice as large. This reflects several factors of 2 
and 4: The factor of (1/2) in the constant C; the factor of 4 in Eq. pl| which expresses that each initial and final state 
appears four times in Eq. the momentum integral for indistinguishable particles extends only over the hemisphere 
whereas distinguishable particles have twice the number of possible final states. 

We used images like Fig. || (for indistinguishable atoms) and Fig. ||a for impurity atoms (see following section) to 
determine the collision cross sections. The cross-section for indistinguishable atoms was found to be 2.1 ± 0.3 times 
larger than the one for impurities in agreement with the expected factor of two [jnj . 

IV. SUPPRESSION OF IMPURITY COLLISIONS 

The key difference between the scattering of light and massive particles (or impurities) is their energy-momentum 
dispersion relation. The dispersion relation for impurities is Eu = (hk) 2 /2M, whereas for light, Eu = Tike with c 
denoting the speed of light. This difference is responsible for the complete suppression of impurity scattering at low 
velocities. For simplicity, we assume equal mass M for the impurity and condensate atoms. 

Energy-momentum conservation, the S function in Eq. ^, requires E^ — Ek-q = Tilo^ . For impurity particles, the l.h.s. 
is always less than vhq, where v — Tik/M is the initial velocity of the impurities. Thus, collisions with the condensate 
are only possible, if this maximum energy transfer is sufficient to excite a quasi-particle, i.e., v > min(cj^/g) = Vl, 



where i>l is the Landau critical velocity |28| for superfluidity below which no dissipation occurs because the density 
of final states vanishes. 

For the excitation spectrum of the condensate (Eq. [li]), the Landau velocity is the Bogoliubov speed of sound 
Vl = c s = y/ p/M. Impurity particles moving below this speed cannot dissipate energy in collisions. In contrast, for 
photons, dEk / d(hk) =c»c s , and scattering is always possible. 

In the perturbative limit (no stimulation by the final occupation), the total rate of scattering Wtot is given by 
integrating Eq. ^ over all possible momentum transfers. In the derivation above for large incident velocities we 
absorbed the angular integral in the density of states since the scattering in the center-of-mass frame is isotropic. 
Now it is necessary to consider the scattering angles explicitly. Furthermore, we stay in the frame where the condensate 
is at rest. Straightforward transformations lead to 



27r|C| 2 A „, ../fck-q hq 2 B \ 

= (N /V){^) 2 j dq ^^S( q )8^^ 



hq 2 B ^ 
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= (N Q /V) a{r,) v, (27) 

where HQ — Mv(l — l/?7 2 ) is the maximum possible momentum transfer, and 77 = v/c s must be larger than 1. 
The collision cross-section is 17(77) = a^F(rf) where en = 47ra 2 . For 77 < 1, F(rf) — and for 77 > 1, F(rf) — 
1 - I/77 4 — log(77 4 )/?7 2 . 

The suppression factor F(rj) is determined by two factors: the phase space restriction due to the 5 function in Eq. |27], 
and additional suppression at low momentum transfers by the structure factor of the condensate. For decreasing 
velocity, the possible scattering angles 9 become restricted to a forward scattering cone (9 < arccos(l/?7)), which 
shrinks to zero solid angle at the Landau critical velocity (Fig. ||). This reflects that near the Landau velocity, the 
scattered particle has "difficulties" to provide enough energy per momentum to create phonons. The maximum energy 
transfer occurs when the momentum transfer is collinear with the incident velocity, i.e., for small angle scattering 
angles. A graph of the suppression factor as function of impurity velocity is shown in Fig. 




FIG. 6. Momentum transfer in collisions. The polar diagram shows the momentum transfer (in units of the initial momentum) 
vs. scattering angle 6 for different values of 77 = v/c s (10, 3, 1.8, 1.5, 1.3, 1.1). As the impurity velocity approaches the speed 
of sound (77 — > 1), the scattering cone shrinks to zero solid angle. 
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FIG. 7. Suppression of collisions. Shown is the suppression factor as a function of the impurity velocity (normalized by the 
speed of sound c 8 , solid line). The dotted line represents the suppression due to phase-space restriction alone (i.e., setting the 
structure factor S(q) = 1). 



Experimentally, the Landau critical velocity can usually only be observed by moving microscopic particles through 
the supcrfluid which do not create a macroscopic flow pattern. Studies of superfluidity with microscopic objects were 
pursued in liquid 4 He by dragging negative ions through pressurized 4 He |2^]3(J, and by scattering 3 He atoms off 
superfluid 4 He droplets pl| . 

To study the effects of impurities interacting with the condensate, we created microscopic impurity atoms using 
a stimulated Raman process which transferred a small fraction of the condensate atoms in the \F — l,mj? = — 1) 
hyperfine state into an untrapped hyperfine state \F — 1, nip — 0) with a well-defined initial velocity [p2| . The initial 
velocity could be adjusted between zero and two single-photon recoil velocities by varying the angle between the two 
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Raman beams. As these impurities traversed the condensate, they collided with the stationary condensate, resulting 
in a redistribution of the impurity momenta which was detected by a time-of- flight analysis |Ll[| (Fig. 0) . 
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FIG. 8. Observation of elastic collisions between the condensate and impurity atoms, (a) Impurities traveling at 6 cm/s 
along the radial axis (to the left in images) were scattered into an s-wave halo. Absorption image after 50 ms of time-of-flight 
shows the velocity distribution after collisions between the condensate (bottom) and the outcoupled tuf = atoms (top), 
spatially separated by a Stern-Gerlach type magnetic field gradient. The collisional products are distributed over a sphere in 
momentum space. The image is 4.5 x 7.2 mm. (b) Similar image as (a) shows the collisional products (arrow) for impurity atoms 
(top) traveling at 7 mm/s along the condensate axis (upward in image). For this image, v g /c s = 2.7 (see text). Collisions are 
visible below the unscattered impurities, (c) Similar image as (b) with v g /c 3 — 1.6. Collisions are suppressed. The momentum 
distribution of the outcoupled atoms was distorted by mean-field repulsion. The images are 2.0 x 4.0 mm. Figure is taken 
from Ref. (jl). £ 
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FIG. 9. Onset of superfluid suppression of collisions. Shown is the observed collisional density normalized to the predicted 
one in the limit of high velocities as a function of r) = v g /c s , which is a measure of the impurity velocity in units of the 
condensate's speed of sound. The error bars represent the statistical uncertainty. Data is taken from Ref. 
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To probe for the suppression of collisions, the impurity velocity has to be varied around the speed of sound. For that, 
we produced impurity atoms at low velocities (7 mm/s) and varied the speed of sound by changing the condensate 
density. The small axial velocity imparted by Raman scattering allowed us to identify products of elastic collisions 
in time-of-flight images (Fig. ||b, c) since collisions with the stationary condensate redistributed the impurity atoms 
toward lower axial velocities. However, the impurity velocity was predominantly determined by the gravitational 
acceleration g, which imparted an average velocity of v g = \[gT z where l z is the Thomas- Fermi diameter of the 
condensate in the z-direction. Thus, the effect of superfluidity on impurity scattering depends primarily on the 
parameter rj — v g /c s which is the ratio of the typical impurity velocity v g to the speed of sound c s at the center of 
the condensate. 

A time-of-flight analysis of impurity scattering for the case of a low-density condensate (small c s ) and large con- 
densate radius (large v g ) is shown in Fig. IJd. The effect of collisions is clearly visible with about 20% of the atoms 
scattered to lower axial velocities (below the unscattered impurities in the image). In contrast, in the case of tight 
confinement, the condensate density is higher (larger c s ) and its radius is smaller (smaller v g ), and the collision 
probability is greatly suppressed due to superfluidity (Figs. |p and [)]). 

The bosonic stimulation factor in Eq. ^ becomes relevant if the final states are populated, either by scattering 
or thermally. We observed that the fraction of collided atoms increased with the number of outcoupled impurities 
(Fig. [l(]). For a large outcoupled fraction, population nk~ q and N q is built up in the final states and stimulates 
further scattering. This collisional amplification is not directional, and is similar to the recently observed optical 
omnidirectional superfluoresence [p3[. 
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FIG. 10. Collective amplified elastic scattering in a Bose-Einstein condensate. Shown is the fraction of collided atoms vs. 
the number of outcoupled atoms. For this data, v g /c a =4.9 and the chemical potential was 1.8 kHz. Figure is taken from Ref. 
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Gain of momentum and thus transfer of energy from the condensate to the impurity atoms is impossible at zero 
temperature, but may happen at finite temperature due to the presence of thermal excitations (the N_ q term in 
Eq. H)El Thus finite temperature enhances the elastic cross-section by two effects: Absorption of quasi-particles 
(anti-Stokes process) and stimulation of momentum transfer by the final state population (Stokes process). 

Fig. [ll] shows the dramatic variation of the elastic scattering cross-section with temperature. However, the finite 
temperature did not affect our data in a major way: Due to gravitational acceleration we couldn't probe the velocity 
regime well below the Landau critical velocity where only thermally assisted collisions are possible. Furthermore, 
when we counted the number of collided atoms we had to use a background subtraction method where we subtracted 
the small signal of the energy gain collisions from the energy loss collisions (see Ref. jll| for details), thus cancelling 
most of the finite-temperature effects. 
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FIG. 11. Temperature dependent cross-section vs. impurity velocity. Shown is the cross-section at zero temperature (solid 
line) and at a finite temperature kT — (i which is typical for our experimental conditions (dash-dotted line). The finite 
temperature cross-section includes collisions involving thermally occupied quasi-particles where the impurities lose or gain 
energy. In the experiment, we measured the number of impurities which lost its energy minus the number which gained energy. 
Thus, the experimental measured cross-sections (Fig. |^) should be compared to ffooii.ioss — o" co n igain (dashed line). 

V. SUPPRESSION OF DISSIPATION FOR A MOVING MACROSCOPIC OBJECT 



So far, we have discussed the suppression and enhancement of microscopic processes (light scattering and impurity 
collisions). The suppression of dissipation is even more dramatic on the macroscopic scale. The flow of liquid 4 He and 
the motion of macroscopic objects through it are frictionless below a critical velocity [[34| . Recently, we have explored 
such frictionless flow in a gaseous BEC Jl^jnjfl , 

The microscopic and macroscopic cases bear many parallels. The onset of scattering or dissipation has two require- 
ments: one needs final states which conserve energy and momentum, and an overlap matrix element which populates 



2 We are grateful to S. Stringari for pointing out the importance of finite-temperature effects. 
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these states. In the case of macroscopic flow, the first requirement leads to a critical velocity for vortex creation and 
the second requirement addresses the nucleation process of vortices. 

The Landau criterion for superfluidity shows that excitations with momentum p and energy E(p) are only possible 
when the relative velocity between the fluid and the walls or a macroscopic object exceeds the Landau critical velocity 
vl which is given by vl — min(i£(p)/p) (see e.g. p|,^4[). A similar criterion applied to vortex formation yields 

£V " M " X 1 \k("] (28) 



MD 



where / V ortox — J P d 3 ^ is the integrated momentum of the vortex ring or line pair, E VOTtcx is its total energy, D is the 
dimension of the container, and £ the core radius of a vortex which in the case of dilute gases is the healing length 
£ = l/y/8irpa. Ref. Q derived Eq. ^8] for vortex rings with a maximum radius D, Ref. looked at pairs of line 
vortices at distance D. Feynman found a similar result for superflow through a channel of diameter D. 

An analogous result is obtained for a Bose condensed system placed under uniform rotation with angular velocity 
Q. A vortex becomes energetically allowed when its energy E' in the rotating frame drops to zero, 



E' = E — VLL = 0. 



(29) 



where E and L are the energy and angular momentum in the laboratory frame. This defines a critical angular velocity 
below which a vortex cannot be sustained due to conservation of angular momentum and energy |37| : 



n c = 



E v 
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(30) 



The critical velocity at the wall of the rotating container, v c = D£l c , agrees with Eq. [28L However, Eqs. |28| and 
|30| only reflect the energy and momentum required to generate vortices, and do not take into account the nucleation 
process. If the scattering particle is macroscopic in size, the coupling is between the ground state and a state containing 
a vortex. Populating such a state requires nucleation of the vortex by the perturbing potential, which usually does not 
occur until higher velocities are reached than those predicted by Eqs. ^8] and [30|. The other option, the formation of the 
vortex by macroscopic quantum tunneling between the two states is an extremely slow process. In recent experiments 
in which a Bose condensate was placed in a rotating potential, the critical angular velocity for the formation of a 
single vortex was observed to be 1.6 times higher than the value given by Eq. |3(] |58|. This discrepancy may be due to 
a nucleation barrier associated with the excitation of surface modes, as some authors have recently suggested p9| , p0| . 

To study frictionless flow in a Bose-condensate, we focused an argon ion laser beam (at 514 nm) onto the condensate, 
which repelled atoms from the focus. The laser beam was scanned back and forth along the axial direction of the 
condensate, creating a moving "hole" that simulated a macroscopic object. Rapid sequence phase-contrast imaging 
allowed us to directly measure the density profile of the superfluid around the moving laser beam. 

For a weakly interacting Bose-condensed gas at density p(r) and chemical potential /i(r), pressure is identical to 
the mean-field energy density P = p(r)p(r)/2 [Q. A drag force arises due to the pressure difference across the moving 
object. The chemical potential is given by p(r,t) = gp(r,t), where g — CV — 4Trh 2 a/M is the strength of two-body 
interactions. The drag force F is given by 



F ~ gSpAp = SpAp/g 



(31) 



where Ap and Ap are the differences in density and chemical potential across the stirring object, and S the surface 
area the macroscopic object presents to the condensate. 
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FIG. 12. Pressure difference across a laser beam moving through a condensate. On the left side in situ phase contrast images 
of the condensate are shown, strobed at each stirring half period: beam at rest (top); beam moving to the left (middle) and to 
the right (bottom). The profiles on the right are horizontal cuts through the center of the images. The stirring velocity and 
the maximum sound velocity were 3.0 mm/s and 6.5 mm/s respectively. Figure is taken from Ref. Jl3| ]. 

If the laser beam is stationary, or moves slowly enough to preserve the superfluid state of the condensate, there will 
be no gradient in the chemical potential across the laser focus, and therefore zero force according to Eq. ||l]. A drag 
force between the moving beam and the condensate is indicated by an instantaneous density distribution p(r, t) that 
is distorted asymmetrically with respect to the laser beam. Fig. |lj shows phase contrast images strobed at half the 
stirring period (where the laser beam is in the center of the condensate). A bow wave and stern wave form in front 
of and behind the moving laser beam, respectively. We define the asymmetry A as the relative difference between the 
peak column densities in front (pf) and behind (p b ) the laser beam A = 2(pf — p b )/{pf + p b ). The asymmetry A is 
proportional to the drag force F. 
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FIG. 13. Density dependence of the critical velocity. The onset of the drag force is shown for two different condensate 
densities, corresponding to maximum sound velocities of 4.8 mm/s (•, left axis) and 7.0 mm/s (x, right axis). The stirring 
amplitudes are 29 pm and 58 pm, respectively. The two vertical axes are offset for clarity. The bars represent statistical errors. 
Figure is taken from Ref. Ea |. 

In Fig. [l3| we show measurements of the asymmetry for two maximum densities po of 9 x 10 13 and 1.9 x 10 14 cm~ 3 . 
In each data set there is a threshold velocity v c below which the drag force is negligible, and this threshold increases 
at higher density. Its value is close to 0.1 c s for both data sets, where c s is the sound velocity. Above this critical 
velocity, the drag force increases monotonically, with a larger slope at low density. 

One can compare measurements of the asymmetry (proportional to the drag force F) with the rate of energy 
transferred to the condensate, F • v, using the calorimetric technique introduced in p"2| . For this, the condensate was 
stirred for times between 100 ms and 8 s, in order to produce approximately the same final temperature. After the 
stirring beam was shut off, the cloud was allowed to equilibrate for 100 ms. Using ballistic expansion and absorption 
imaging, we determined the thermal fraction and thus the temperature and total energy. 
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FIG. 14. Calorimetry of a condensate. The energy transfer rate during stirring (•, left axis) was obtained from temperature 
measurements. The error bars reflect shot-to-shot variations in the temperature. The results are compared to the energy 
transfer rate vA(v) obtained from asymmetry measurements of the flow field during the stirring (x, right axis). Figure is taken 
from Ref. @. 

The calorimetric measurements can be compared with the drag force inferred from the asymmetric density distri- 
bution. Using Eq. Bll the energy transfer rate per atom is written in terms of the asymmetry as 



dE 



8 iiqPqI z D 

vA(v) (32) 



m N 15 N 

where D is the diameter of the laser beam and l z the Thomas Fermi diameter in the radial direction 



Fig. 14 shows that the calorimetric and the drag force measurements are in remarkable agreement over the entire 
velocity range up to a single scale factor for vA(v), demonstrating the consistency between the two methods. For the 
parameters of our experiment (D ~ 10/im, po — 1.3 x 10 14 cm~ 3 , l z = 66/im, N = 1.8 • 10 7 ) the overall heating rate 
predicted by Eq. ^ is 2.4 times larger than that obtained directly from calorimetry. Possible explanations for this 
difference are the inhomogeneous density profile or a breakdown of the quasi-static approximation embodied in Eqs. 
|l| and ||. 

The observed critical velocity may be related to the formation of vortices. An estimate based on Eq. [2^ for 
typical experimental parameters in sodium (D = 10 /im, peak density po = 1.5 x 10 14 cm~ 3 , a = 2.75 nm) yields 
v c ~ 1.0 mm/s, close to the experimental observations. However, Eq. depends only weakly on the speed of sound, 
through the logarithmic dependence on the healing length £. In contrast, our measurements show an approximate 
proportionality to the sound velocity jl3| , suggesting that vortex nucleation determines the onset of dissipation. 

Time-dependent simulations of the Gross-Pitaevskii equation show the formation of vortex line pairs above a critical 
velocity which is close to the observed value 0] . Several authors have emphasized the role of locally supersonic flow 
around the laser beam in the nucleation of vortices |^l|-^3|. In one theoretical model [0, the vortices are emitted 
periodically at a rate that increases with velocity, and reduce the pressure gradient across the object. The predicted 
heating rate [^DjiJ] is in rough agreement with the experimental results. Moreover, this model also predicts that the 
slope of the asymmetry vs. velocity should increase at lower density, in accord with our observations (Fig. |l3| ). 

The calorimetric measurements were extended to purely thermal clouds |j45"f . By accounting for the different 
geometries of the clouds, we could infer the effective energy transfer to an atom in a collision with the stirrer. For 
velocities well above the critical velocity we found almost equal energy transfer. This result fits well into the picture 
that the special properties of a condensate only show up in processes where little energy is transferred to the atoms. 



VI. FOUR- WAVE MIXING OF LIGHT AND ATOMS 



Our discussion on amplification of light and atoms in a BEC in the next section will reveal novel aspects of the 
coupling between light and atoms. In preparation for this, we want to present some general aspects of the scattering 
Hamiltonian in Eq. [TJ applied to the scattering of photons from the condensate (the operators c\ , Ck are now creation 
and annihilation operators for photons). The coupling can be regarded as four- wave mixing of two atomic fields and 
two electromagnetic fields. 

When an atom is illuminated by two strong laser beams in a A configuration similar to Fig. [IJ the electronically 
excited state can be adiabatically eliminated for sufficiently large detuning A. The coupling matrix element between 
the two atomic ground states is Tifl R /2, where flu is the two-photon Rabi frequency. This Rabi frequency can 
be expressed by the (complex) electric field strength Eip of the two laser beams: flu — flifl2 cos 0/2A, where 
fli,2 = Ei^d/h are the Rabi frequencies of the individual laser beams with the atomic dipole matrix element d. <fi 
is the angle between the axes of polarization of the two laser beams. The two-photon Rabi frequency flu can be 
rewritten as 

fl R = d 2 E 1 E 2 cos (f>/2h 2 A. (33) 

Using the Hamiltonian in Eq. |l|, the coupling matrix element squared is \C\ 2 n\n2 where are the photon numbers 
in the two beams. The two-photon Rabi frequency fl R is given by 

(hfl R /2) 2 = \C\ 2 n x n 2 . (34) 
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Expressing the number n± of photons with angular frequency too in a volume V by the complex electric field strength 
£1 



m = e |£i| 2 U/2^ , (35) 
and comparing to Eq. |3^ we obtain for the coupling constant C between the two modes 

w d 2 cos </> 

c = ^vA- (36) 

Another simple limit of the Hamiltonian in Eq. [I] is the situation when atoms are illuminated by a single laser beam 
in mode 1. Then the diagonal term Ccjajciao gives rise to the AC Stark shift. With the photon number ni, the AC 
Stark shift AE of an atom is AE = Gn\. Using Eqs. [35] and ^ with cos0 = 1 one obtains the well known result, 
AE = d 2 E\ /Ah A. 

Since the four- wave mixing Hamiltonian (Eq. |l]) applies to interactions of a condensate both with light and atoms, 
we can draw analogies. The AC Stark shift corresponds to the mean field interaction between impurity atoms and the 



condensate. In Sec. |Hj we saw how the scattering of atoms into empty modes gave rise to the usual elastic collision 
rate. Similarly, the scattering of photons into empty modes by the four-wave mixing Hamiltonian results in Rayleigh 
scattering which we want to discuss now in more detail. 

Rayleigh scattering is described by Eq. 0. In the limit of weak scattering and for non- interacting atoms (S(q) = 1) 
one obtains 

W+/N = T \C\ 2 ni 8(E k - E k _ q - hio° q ). (37) 

The scattering rate per solid angle d7 sca tt/df2 is given by integrating this equation over all final states using the 
density of states per energy interval and solid angle 

dp VUJ Z (38) 



dEdfl {2irc) 3 h 

Using Eqs. |35|, |36| one obtains 

d7scatt _ 27T 2m dp _ uj^Efd 4 cos 



dfl h 1 1 L dEd£l 32TT 2 h i c 3 e A 2 
With the expression for the natural linewidth T 



tIpiv^= " 3 , (39) 



_fu4_ 



(40) 



this simplifies to 



d7 S catt 3cos 2 </> 3sin 2 6» 

1 = — z ^' v 41 ) 



dQ 8tt 4A 2 8tt 

where we have defined the Rayleigh rate R 

R =^- ( 42 ) 

9 denotes the angle between the linear polarization of the incident light and the direction of the scattered light. For 
each scattering angle 6 there are two polarizations of the scattered light. One is orthogonal to the incident polarization, 
thus cos0 = in Eq. [36], and it doesn't contribute. The density of states in Eq. [38] was therefore defined for only 
one polarization. The other polarization is in the plane of the incident polarization and the scattering direction, thus 



cos^ = sin6>. R was defined in Eq. 42 in such a way that integration of Eq. |41| over the whole solid angle gives 

7scatt = R- 

Below, we will need another useful expression for the two-photon Rabi frequency We use Eq. |34|, assume that 
one laser beam is a weak probe beam with photon number n p , express the photon number of the other (strong) beam 
by its Raleigh scattering rate and obtain 
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n 2 R = 6nRX^cn p /V. (43) 

Four-wave mixing between atoms and light has two important limiting cases. In one case, atoms are diffracted by 
a standing wave of light, i.e., the atoms move in the AC Stark shift potential of two interfering light fields. In the 
other case, light is diffracted by a matter wave grating, a density modulation formed by two interfering matter waves. 

In the first case, the AC Stark shift potential is proportional to Jn^ni , where n k ,i are the number of photons in the 
two laser beams. The diffraction efficiency is proportional to the square of the potential, and therefore, for N q atoms 
in the incident mode, the scattering rate is proportional to HkniN q . In the latter case, the density modulation caused 
by the interference between N k and Ni atoms is proportional to y/NkNi. The scattering rate for light with n q photons 
is then proportional to NkNin q . Applying this to a condensate with No atoms illuminated with a (strong) laser beam 
with nfc photons, we see from Eq. ^| that for N q < rik- q one has Bragg scattering, the scattering of atoms from a 
standing wave of light. This is usually realized by illuminating the atoms with two laser beams. For N q > n k - q the 
physical picture is the diffraction of light by an atomic density modulation. It is this regime, which we have exploited 



for the amplification of atoms in a Bose-Einstein condensate (see Sect. VIII for a further discussion). 

Generally, the optical (or atomic) gratings are moving. As a result, the diffracted atoms (or photons) have an 
energy different from that of the incident particles. However, one can always transform to a moving frame where the 
grating is stationary and there is no energy transfer in the scattering process. The frequency shift due to the moving 
grating can therefore be regarded as the Doppler shift related to the Gallilean transformation between the two frames. 

VII. SUPERRADIANCE AND MATTER WAVE AMPLIFICATION 

Spontaneous light scattering can be stimulated when the atomic recoil state is already populated (the N q term in 
Eq. 0). We have explored this process in our studies of superradiance fljij , phase-coherent atom amplification [ JL7| , 
and optical amplification Jl8| . 

In all these experiments, the condensate was illuminated with a laser beam (mode fc, also called the "dressing 
beam"). A condensate atom scatters a photon from the laser beam into another mode and receives the corresponding 
recoil momentum and energy. Injection of atoms turns this spontaneous process into a stimulated process and realizes 
an amplifier for atoms. The injected atoms interfere with the condensate at rest and form a matter wave grating 
which diffracts the dressing light. The diffraction transfers recoil momentum and energy to the atoms, which results 
in a growth of the grating and therefore the number of atoms in the recoil mode — this is the intuitive picture for atom 
gain. If no atoms are injected, the whole process may start from spontaneous scattering as superradiance. We will 



discuss in Sec. VIII that the build-up of the matter wave grating can be induced also by a probe light beam resulting 
in optical amplification. 

Eq. [7] describes the scattering rate out of a condensate into a recoil state with population N q . In the limit of an 
empty mode for the scattered light (nk- q = 0), it reduces to 

W+ = ^-\C\ 2 N n k (N q + l)6(E k ~ E k _ q - Huf). (44) 

For the high momentum transfers considered here (on the order of the photon recoil momentum), S(q) = 1. Each 
scattering event which transfers momentum ?iq to the condensate, generates a recoiling atom in mode q. The final 
states of the photon which are associated with a momentum transfer hq form a continuum. Integrating W + over all 
such final states gives the growth rate N q for the recoiling atoms. Using Eq. [IT] we obtain 

N q = G q {N q + l)-T 2 , q N q (45) 

with the gain coefficient 

G q = RN S ^Q q . (46) 

The solid angle £l q reflects the number of photon modes which are excited together with a quasiparticle with 
momentum ?iq. Nq is the number of atoms in the condensate at rest and 8 q is the angle between the polarization of 
the dressing beam and the direction of photon emission. In addition, a loss term T2, q was included which describes 
the decoherence rate of the matter-wave grating and determines the threshold for exponential growth. It represents 
the linewidth of the two-photon process which generates recoil atoms in mode q. The total scattering rate Wt Q t is 
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(47) 



It is important to realize how the solid angle is divided into a sum over modes q. When iV g <l, the sum in Eq. H 
is simply an integral over the whole solid angle and the result (without the loss term) is the Rayleigh scattering rate: 
Wtot = = RNo. However, when the build-up of population N q becomes important, the division of the solid angle 
into "coherent pieces" becomes essential. The light, emitted into the direction k— q, and the quasi-particles are created 
in the finite volume of the condensate. Therefore, their momentum is only defined to within h over the dimension of 
the condensate. Each "mode" represents a solid angle corresponding to that uncertainty (the longitudinal momentum 
of the photon is determined by energy conservation, i.e., the S function in Eq. Q). Therefore, each scattering mode 
spans a diffraction limited solid angle which is approximately fl q ss A 2 / A q where A is the wavelength of the scattered 
light and A q is the cross-section of the condensate perpendicular to the axis of light emission. 

Our Les Houches notes discuss a semiclassical derivation of the gain mechanism |Tc|] . In this picture, one has 
N q atoms in the recoil mode interfering with the condensate at rest resulting in a modulated density. When this 
atomic distribution is illuminated with the dressing light, all the atoms can be regarded as driven oscillators. Their 
radiation interferes constructively in the direction 6 q in which the density grating diffracts the dressing light. The 
phase-matching condition is fulfilled for this direction and a solid angle fl q around it. This angle can be rigorously 
obtained from the usual phase-matching integral for superradiance in extended samples Eq] 



a, 



dfi(k) 



p(r) exp(i(ki — k) • r)dr 



(48) 



where k^ is the wave vector of the incident light, |k| = |k»|, and p(r) is the normalized atomic density distribution 
(fp(r)dr = l). 

The key results of the discussion above are the gain equations, Eqs. |4g and [46]. For a condensate of cross-section 
A q and length l q , the gain is proportional to Noil q w po\ 2 l q which is proportional to the resonant optical density of 
the condensate along the direction of the scattered light. Therefore, for an anisotropic condensate, the gain is largest 
when the light is emitted along its longest axis (the "end-fire mode" fl47|l ). When the intensity of the dressing light is 
above threshold (G q > I^q in Eq. |i|), a condensate will emit highly directional beams of light and atoms along the 
direction for which the gain is highest. Our observation of this phenomenon has been described in Refs. |Tc|Jl^ ]. 




FIG. 15. Input-output characteristic of the matter-wave amplifier, (a-c) Typical time-of-flight absorption images demon- 
strating matter wave amplification. The output of the seeded amplifier (c) is clearly visible, whereas no recoiling atoms are 
discernible in the case without amplification (a) or amplification without the input (b). The size of the images is 2.8 mm x 
2.3 mm. (d) Output of the amplifier as a function of the number of atoms at the input. A straight line fit shows a number 
gain of 30. Reprinted by permission from Nature, Ref. flT] , copyright 1999 Macmillan Magazines Ltd. 



The amplification of atoms is conceptionally even simpler: Atoms in a certain recoil mode q are injected into the 
condensate and are amplified with a gain coefficient G q . We only summarize the main experimental results — for a 
full account see Refs. ]fp| , p7t . 
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FIG. 16. Experimental scheme for observing phase coherent matter wave amplification. A small-amplitude matter wave 
was split off the condensate by applying a pulse of two off-resonant laser beams (Bragg pulse). This input matter wave was 
amplified by passing it through the condensate pumped by a laser beam. The coherence of the amplified wave was verified by 
observing its interference with a reference matter wave, which was produced by applying a second (reference) Bragg pulse to 
the condensate. The total number of atoms in the recoil mode showed constructive and destructive interference between the 
amplified input and the reference matter wave as the phase of the reference wave was scanned. Reprinted by permission from 
Nature, Ref. |l7|, copyright 1999 Macmillan Magazines Ltd. 



Input matter waves with momentum hq were produced by exposing the condensate to two laser beams. Their 
difference frequency was tuned to the resonance for Bragg scattering (see Sec. [fl]). The stimulated redistribution of 
photons among the two beams transferred recoil momentum to the atoms. The fraction of atoms in the recoil state 
was controlled by the intensity and duration of the Bragg pulse. While these input atoms were still in the condensate 
volume, they were amplified when the condensate was exposed to the dressing beam. Fig. [l5] shows the input-output 
characteristics of the amplifier. The gain was controlled by the intensity of the pump pulse (see Eq. |46| ) and typically 
varied between 10 and 100. Fig. |l5| d shows the observed linear relationship between the atom numbers in the input 
and the amplified output with a number gain of 30. 

This atom amplifier is a narrow band amplifier. It only amplifies input momentum states which can be populated 
by condensate atoms by scattering a photon of the dressing beam. The possible input states lie on a sphere in 
momentum space. Its center is displaced from the origin (the momentum of the condensate at rest) by the momentum 
Jik of the dressing beam. The thickness of the momentum sphere is determined by the momentum uncertainty of the 
condensate, which is Ti over its size, as was directly measured using Bragg spectroscopy [^]. In our experiment, the 
input momentum was automatically matched to the amplifier's narrow bandwidth since the input beam was created 
by Bragg scattering, and one of the Bragg beams was identical to the dressing beam. 

The Hamiltonian in Eq. [l] provides coherent coupling between the light and atoms. Therefore, the atom amplification 
should be phase-coherent. This was experimentally verified with an interferometric technique. For this, a reference 
matter wave was split off the condensate in the same way as the first (input) wave (see Fig. |l6|). The phase of 
the reference matter wave was scanned by shifting the phase of the radio-frequency signal that drove the acousto- 
optic modulator generating the axial Bragg beam. We then observed the interference between the reference and the 
amplified matter waves by measuring the number of atoms in the recoil mode. 



VIII. AMPLIFICATION OF LIGHT IN A DRESSED CONDENSATE 

A dressed condensate (Fig. [n]), a condensate illuminated by laser light, was used as an atom amplifier. Now we 
develop this picture further. Rayleigh scattering produces scattered photons and recoiling atoms. In the dressed atom 
picture, this is described as the decay of the dressed condensate into a photon and recoiling atom, or in other words, 
the dressed condensate can spontaneously emit pairs of photons and atoms. The amplification of atoms discussed in 
the previous sections solely focuses on the recoiling atoms "emitted" by the dressed condensate. Although recoiling 
atoms and scattered photons are emitted in pairs, the photons leave the condensate almost instantaneously and there 
is no significant population build-up. Formally, as discussed in the Les Houches notes, one can adiabatically eliminate 
the light field from coupled equations and obtain the gain equation for the matter waves. 

On the other hand, the dressed condensate should act also as an amplifier for light. An input optical field should 
stimulate Rayleigh scattering processes which results in photons scattered into the input mode. Our recent experiments 
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on optical amplification in a BEC fig ] required a more general description of the interplay between optical and matter 
wave amplification. 

On the following pages, we present a general discussion of four-wave mixing of light and atoms. We first start 
simply with the gain cross-section for the light and a complex index of refraction. The recoiling atoms enter the 
picture in two stages, first within the framework of Heisenberg equations in the undepleted-pump approximation, and 
then using optical Bloch equations. 



A. Cross-section for optical gain and slow light 



The physical picture behind the optical gain of the dressed condensate is as follows: if a very weak probe beam is 
injected into the dressed condensate, it acts together with the dressing beam as a pair of Bragg beams and creates 
recoiling atoms. This process transfers photons from the dressing beam into the probe beam. At higher gain, the 
recoiling atoms become significant. They move out of the condensate (or decohere) on a time scale T^" 1 which is the 
inverse of the linewidth of the Bragg transition. In steady state, the number of recoiling atoms N q in the volume of 
the condensate is proportional to the intensity of the probe light. Those recoiling atoms interfere with the condensate 
at rest and form a diffraction grating which diffracts the dressing beam into the path of the probe light resulting in 
amplification of the probe light (Fig. 17). 
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FIG. 17. Amplification of light and atoms by off-resonant light scattering, (a) The fundamental process is the absorption 
of a photon from the "dressing" beam by an atom in the condensate (state |1)), which is transferred to a recoil state (state 
|2)) by emitting a photon into the probe field. The intensity in the probe light field was monitored by a photomultiplier. (b) 
The two-photon Raman-type transition between two motional states (|1), |2)) gives rise to a narrow resonance, (c) The dressed 
condensate is the upper state (|1')) of a two-level system, and decays to the lower state (recoil state of atoms, |2)) by emitting 
a photon. Figure is taken from Ref. [[L8| 



An expression for the gain can be derived in analogy to a fully inverted two-level system with dipole coupling which 
would have a gain cross-section of 67rA 2 for radiation with wavelength A(= 2nX). For the Raman-type system in Fig. 
|l7jb, the gain is reduced by the excited state fraction, R/T (where R is the Rayleigh scattering rate for the dressing 
beam and T is the linewidth of the single-photon atomic resonance) and increased by r/T2, the ratio of the lincwidths 
of the single-photon and two-photon Bragg resonances. Thus the expected cross-section for gain is 



t 2 



(49) 



The lineshape of the optical gain is that of the two-photon Bragg resonance. Due to the long coherence time of 
a condensate, it has a very narrow linewidth. Such a narrow band gain is accompanied by a slow group velocity of 
light. This can be described by a complex index of refraction n(u) = n\{u>) + in2(w). 

For a Lorentzian resonance curve with FWHM of T2 for the gain, the complex index of refraction is 



1 + s 2 1 + s 2 



,9'i 



(50) 



where is the background index of refraction, g' denotes the strength of the resonance, and <5 = (lo — loq) / '(1 2 / '2) is 
the normalized detuning from the resonance at loq. The imaginary part of Eq. ^ has the usual Lorentzian lineshape. 

The real and imaginary parts of the index of refraction are connected by Kramers-Kronig relations. For the special 
case above of a Lorentzian lineshape the gain and dispersion at resonance are connected by 
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UJ — UJq 

A steep slope of the (real part of) the index of refraction gives rise to a slow group velocity of light 

c 



(51) 



(52) 



Eqs. [5I] and [52] imply a simple relationship between the gain and delay time for an optical pulse. The amplitude 
of an optical pulse at frequency ujq which propagates through a medium with the index of refraction n of length I is 
amplified by a factor 



g = cx-p(-n 2 Ld Q l/c). 

When the first term in the denominator of Eq. ^ is dominant, the delay time 773 of the pulse is 

/ too dni 2 In g 
t d = — ~ I— -7— = -pr— ■ 
Vg c dw 1 2 



(53) 



(54) 



This equation provides a simple relationship between a narrow band gain and pulse delay pS| ]. A non-inverted 
absorptive two- level system gives rise to "super luminal" pulse propagation p9| . 

For the experimental study of the optical gain, a condensate was illuminated ( "dressed" ) with a single off-resonant 
laser beam and probed with another laser beam, which was red-detuned by 91 kHz to satisfy the Bragg resonance 
condition. Both the dressing beam and the probe beam were in the plane perpendicular to the long axis of the 
condensate, and intersected at an angle of 135 degrees. The probe beam, which propagated parallel to the axis of 
imaging, was much larger than the condensate size. In order to block all the light that did not pass through the 
condensate, a slit was placed at an intermediate imaging plane. The light transmitted by the slit was recorded with 
a photomultiplier. The polarization of each beam was set parallel to the long axis of the condensate to suppress 
superradiance to other recoil modes jljj . 

Fig. [ll] shows that light pulses were delayed by about 20 fis across the 20 fim wide condensate corresponding to a 
group velocity of 1 m/s. This is one order of magnitude slower than any value reported previously (see Ref. |5(J and 
references therein). Fig. |l8| b presents the experimental verification of the relationship between gain and delay time 
(Eq. g). 
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FIG. 18. Pulse delay due to light amplification, (a) About 20 ps delay was observed when a Gaussian pulse of about 140 fis 
width and 0.11 mW/cm 2 peak intensity was sent through the dressed condensate (bottom trace). The top trace is a reference 
taken without the dressed condensate. Solid curves are Gaussian fits to guide the eyes, (b) The observed delay was proportional 
to the logarithm of the observed gain. Figure is taken from Ref. . 



B. Relation between optical gain and atomic gain 

Both the optical gain g (Eqs. |4^ and |5^) and the matter wave G gain (Eq. [l6|) have the same origin, stimulated 
Rayleigh scattering. Therefore, the two gain coefficients should be related. The expression for G (Eq. ^) involves a 
solid angle factor £l q which is proportional to A 2 /A. To be consistent with the optical Bloch equation (to be discussed 
below) we use now 2A 2 /A for fl q and obtain 

G=^^A 2 = (p Kain/2)r 2 . (55) 
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The gain g for the amplitude of the optical field is 

G 

g = exp(p ^ gain /2) = exp(G/r 2 ) « 1 + — (56) 

1 2 

with the last equation being an approximation for small gain. However, Eq. [56] cannot be universally valid. When the 
gain G is above the threshold for superradiance, G > T 2 (Eq. [l5|) the optical gain should diverge: a single recoiling 
atom created by the probe light and dressing light is exponentially amplified and creates a huge matter wave grating 
which will diffract the dressing light into the probe light path, thus amplifying the probe light by a divergent factor. 
Indeed, as we will derive below, Eq. |5(] is only valid at small values of G/T 2 - 

If we can neglect the depletion of the dressed condensate and the dressing laser beam we can simplify the interaction 
Hamiltonian in Eq. |l| to 

H' = C'{tfc i +ac). (57) 

Here, d(c) indicates the atomic (light) field to be amplified, but the following derivation is completely symmet- 
ric between the two fields. This Hamiltonian is a standard down-conversion Hamiltonian. Here it describes the 
down-conversion of the dressed condensate into photons and recoiling atoms. Considering only two modes neglects 
propagation effects in the amplification. The coefficient C defines the time constant of the amplification and is 
proportional to the amplitude of the dressing beam and the square root of the number of atoms in the condensate. 
The Heisenberg equations of motions are 

ik = [a,H'] = C'c\ (58) 
ic't = [c\H'] = -C'a. (59) 

This leads to exponential growth of a and c (proportional to exp(C'i)). However, to describe the physical situation 
in the experiments, one has to allow for damping by introducing r a (r c ) as phenomenological damping time constants 
for a(c). We also include source terms (input fields) and approximate the operators by c numbers: 

a = — ^-(a - 00) - iC'c*, 

6* = iC'a- J ±(c*-4). (60) 

The solutions show relaxation (C 2 < r a r c /4) or exponential growth (C 2 > r a r c /4) depending on the strength of 
the coupling relative to the damping rates. The "gain" below the threshold can be defined as a(t — > oo)/ao for atoms 
(assuming cq = 0) and as c(t — > oo)/co for light (assuming ao = 0), yielding 

_ a(t -> 00) _ c(t -> 00) r a r e /4 
9 ~ a ~ Co ~ (r a r c /4) - c 2 ' [b ' 

The fact that the two gain coefficients are equal is a general property of parametric amplification where two kinds of 
particles are produced in pairs. 

In the limiting case that one field is strongly damped (e.g. that light quickly escapes from the system, T c ^> r o ), 
one can adiabatically eliminate this field from the coupled equation (assuming no photon input (cq = 0)) 

c* = (62) 

and obtain a single gain equation for a. The gain equation for the atom field is 

T 2C" 2 

a = — ^-(a- a ) + -—a. (63) 
z 1 c 

In the absence of damping, the atom number would increase exponentially with a rate constant 4C" 2 /r c which 
we therefore identify with the atom gain rate coefficient G in Eq. [5^. This can be shown explicitly using C' 2 — 
\C\ 2 N Q nk/h 2 and setting the mode volume V = AL, where A is the cross section of the condensate. The axial length 
L could be the condensate length I, but will cancel out. Eqs. [34| and ^ yield 4C" 2 /r c = 2Gc/LT c which equals G 
when we set the decay rate T c /2 equal to the photon transit time c/L. We can then rewrite the gain calculated above 
as 
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r r 

a = — = — (64) 

9 T a - 4C' 2 /r c T a -G ( ' 

For the dressed condensate, we identify T a with r 2 . As expected, at the threshold to superradiance (G = r 2 ), the 
(steady-state) gain for both light and matter waves diverges. The gain can be rewritten as 

9 = r 2 /(r 2 - G) = 1 + G/(r 2 - G). (65) 

In the low gain limit, this yields the same result as Eq. |5(| The comparison with Eq. [56] shows that the effect of the 
coupled equations is to replace the two-photon linewidth r 2 in Eq. [56] by the dynamic coherence decay rate T 2 — G. 
Since propagation effects have been excluded, we can't expect to obtain the exponential factor in Eq. [56], but rather 
the linearized form. The expansion 

g = l + (G/r 2 ) + (G/r 2 ) 2 + ... (66) 
describes the transition from (linear) single-atom gain to (nonlinear) collective gain. 



C. Optical Bloch equations 

The discussion in the previous two sections assumed that the condensate is undepleted — i.e., we calculated properties 
of a condensate with all the atoms in the initial dressed state. However, the presence of the dressing light and the 
probe light depletes the condensate. Furthermore, the calculated amplification coefficients are only valid in a quasi- 
steady state regime which is usually preceded by transient behavior. A correct interpretation of the experimental 
results required a more complete description of the dynamics of the system which will be developed in this section 
using optical Bloch equations. 

We proceed in two steps. In the limit of weak optical gain (or strong probe laser intensity), we will use the ordinary 
optical Bloch equations where the laser fields are treated as constant. Later we will introduce an additional equation 
for the dynamics of the probe light. The condensate at rest (|1)) and the atoms in the recoil state (|2)) are treated as 
a two-level system coupled by the four-wave mixing Hamiltonian which gives rise to a two-photon Rabi frequency £l R 
(Eq. |3^). The coherence between those two states decays at a rate T 2 /2. Assuming constant fl R , the optical Bloch 
equations at resonance take the following simple form 

r 2 

v = — — v-Q R w (67) 
w = n R v (68) 

where v = 2 Im(/7i 2 ) represents the amplitude of the matter wave grating (/ay is the atomic density matrix) and 
w = P22 — Pn is the population difference between the two states pl[ . 
The eigenvalues of the matrix 



-r 2 /2 -n R 
n R o 



(69) 



are A± = — T 2 /4 ± a/ (r 2 /4) 2 — £l R . In the limits of large and small laser intensities one obtains 



This means that at high intensities the system exhibits damped oscillations — Rabi oscillations between the two levels. 
At low intensities, there is relaxation in two steps: The coherence is damped with a rate of T 2 /2, followed by depletion 
of atoms in the condensate, which happens at a rate of 2£!|,/r 2 . It is in this temporal window (2/r 2 < t < r 2 /251|.) 
that the perturbative treatment with the complex index of refraction applies. For longer times, the condensate 
becomes depleted and the assumption that most of the atoms are in the initial dressed state is no longer valid. 

The optical Bloch equations can be analytically solved for a step function input. With the initial condition that at 
time t = all the atoms are in the condensate at rest (w(t = 0) = —I, v(t = 0) = 0) one obtains 
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v(t) 



,/fi|-(r 2 /4) ; 



t) sin 



(r 2 /4) 2 t) 



exp(-^t) sinh v/( r 2/4) 2 - fi£ * 



(71) 



simplifying in the limit of small probe laser intensity (CIr <C 1^/4) to 

2Cl R 



v(t) 



r 2 
r 2 



-exp(-^t) 
-exp(-yt) 



exp( 



r 2 ; 



t < r 2 /n|. 



(72) 
(73) 



By reducing the probe power, the Rabi oscillations slow down and become overdamped and a (quasi-)steady state 
gain is obtained. Inserting Eq. [72| into Eq. |68|one obtains the transition rate Nqw/2 — NqCIrv/2 which is the number 
of photons per unit time emitted by the dressed condensate. To obtain the gain one has to normalize by the input 
photon flux cn p / 1 where n p is the number of photons in the condensate. The amplitude gain is then (assuming small 
gain) 
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1 , N n R i 

1 H ; V. 



(74) 



Using the asymptotic behavior of Eq. |7^ (v(t) ~ 2CIr/T2), Eq. for flu and Eq. 49 one obtains g = 1 + paa m i u l/2 
which agrees with Eqs. |56| and ^5] in the low-intensity limit. Eq. |74| thus has the correct asymptotic limit, but it also 
describes transient behavior when the general solution for v(t) (Eq. |7l|) is used. Theoretical traces based on Eq. |74| 
are directly compared to the experimental results in Fig. |l9[ 

In the experiment, we used long square probe pulses for the probe light (Fig. |l^). When the dressing beam was 
suddenly switched off, a sudden change in the observed probe light intensity was evidence for optical gain. At the 
lowest probe intensity, the depletion of atoms in the condensate was negligible and a clear step at the switch off was 
observed, corresponding to a gain of w 2.8. The initial rise time of « 100 fxs is the coherence time of the dressed 
condensate. At high probe laser power we observed Rabi oscillations in the transmitted probe light. Note that all 
the traces were normalized by the probe beam intensity, and the oscillatory trace at the bottom was obtained at the 
highest probe beam intensity. The oscillations reflect simple two-level Rabi oscillations of atoms between the two 
motional states driven by the two-photon Bragg coupling. 
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FIG. 19. Gain and temporal behavior of light pulses propagating through a dressed condensate, (a) Observed probe pulse 
output from a dressed condensate. The probe light intensities were 5.7mW/cm 2 (bottom), 1.5mW/cm 2 (middle), 0.10mW/cm 2 
(top), while the dressing beam intensity was 5mW/cm 2 , which was just below the threshold for superradiance. The plotted 
signals were normalized by the incident probe intensity and show the gain for the probe light, (b) Calculated probe light 
output for typical experimental parameters. Rabi oscillations develop into steady state gain as the intensity of the probe light 
is reduced. Figure is taken from Ref. Jl^ ]. 

When the probe laser frequency was detuned from the two-photon resonance, the frequency of the Rabi oscillations 
increased (Fig. 20). Optical Bloch equations with detuning [SI predict oscillations at 



flog = J CI 



Au 



(75) 



where Aw is the detuning, in agreement with observations. 
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FIG. 20. Rabi oscillations of atoms observed by monitoring the transmitted probe light intensity, a) Rabi oscillations for 
different detunings of the probe light from the two-photon resonance, b) Observed frequencies of Rabi oscillation versus 
detuning. The solid line is a fit to Eq. c) Amplitude of the oscillations versus detuning. The solid line is a Lorentzian fit 
with a linewidth of about 4 kHz. 



For large optical gain, the Rabi frequency Q R increases during the pulse and the above treatment is no longer valid. 
Therefore, we derive now a second equation which treats the Rabi frequency as a dynamic variable. The population 
transfer to the recoil state (w) results in an increase of the number of the probe beam photons inside the condensate 
volume 

ri p = c(n° - n p )/l + N w/2, (76) 

where I is the length of the condensate with No atoms and crip /I is the input photon flux. Without gain, the steady- 
state number of photons in the condensate volume would be n®. Eq. |76| neglects propagation effects of the light by 
replacing the non-uniform electric field by an average value (such a "mean-field model" is only qualitative, see J52J). 
Replacing the photon number by the Rabi frequency (Eq. |4^) leads to 

2Q R fi R = y(ng - n 2 R ) + ^^^n RV (77) 

where fio is the two-photon Rabi frequency due to the input probe beam and the dressing beam. For small gain, we 
approximate — ^r ~ 2f2fl(f2rj ~~ &r)- This approximation should retain the qualitative features of the coupled 
light-atom system even when the small gain approximation is no longer quantitative. Indeed, we will obtain results 
consistent with our previous treatment (Eq. |65| ) which was not limited to small gain. Using Eq. |5^ for the atom gain 
G we obtain 

ft'ii=|(fio-fiii+yv) (78) 

This equation together with Eqs. ^ and ^ forms a set of coupled equations describing the combined dynamics 
of the atom and light fields. The situation is analogous to the optical laser, where the atomic polarization and the 
electric field inside the cavity are coupled. However, the role of atoms and light is reversed: in the optical laser, the 
cavity lifetime is usually longer than the coherence time of the atomic polarization, whereas in our case the extremely 
long coherence time of the condensate dominates. This would correspond to the bad cavity limit of the optical laser 
which is usually not realized (see J53[ and references therein). 

Assuming rapid relaxation of the light field (ti,R = in Eq. [78]) leads to 

fl R = n + -v. (79) 

Inserting this into Eqs. ^ and ^38] adiabatically eliminates the light field. This treatment is more general than the 
Heisenberg equations above, where we had neglected condensate depletion. To check for consistency, we now assume 
an undepleted condensate (w = —1) and obtain 

v = ^—^v + n . (80) 
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Below the threshold for superradiance, (G < T2), v relaxes with a time constant of 2/(^2 — G) to v — 2f2o/(r2 — G). 
This and Eq. [79 show that the gain g for the probe beam is 



in agreement with Eq. pi 



D. Optical probe of matter wave gain 

The matter wave grating formed inside the condensate is responsible for both atomic and optical gain. We now 
briefly describe experiments where the dynamics of the matter wave grating could be directly observed by monitoring 
the probe light. We first created a matter wave grating with a Bragg pulse and then observed its time evolution by 
monitoring the diffracted dressing beam. The initial seed pulse was 100 us long and transferred about 5% of the 
atoms to the recoil state. 

At lower intensities for which atom amplification was negligible, the grating showed a simple decay (Fig. |l|). At 
higher intensities, collective gain started to compensate the loss, and at intensities above a threshold, net amplification 
was observed. The initial growth rate (Fig. ETJ) followed the linear dependence on the intensity of the dressing beam 
(cx (G — T2)) predicted by Eq. ^ and Refs. |14],Q . The net growth of the matter wave grating was studied previously 
by observing an increase in the number of recoiling atoms in time-of- flight images [ [r7| , whereas Fig. |^ was obtained 
by monitoring the dynamics of amplification in situ by observing light instead of atoms. 

Extrapolating the decay rate in Fig. |l]to zero intensity of the dressing beam gives the decay rate of the matter wave 
grating T2 of (100 /is) -1 , in fair agreement with the linewidth of the Bragg excitation process observed previously [ p0[ . 
This observation of the decay of the matter- wave grating can be regarded as pump-probe spectroscopy of quasi-particles 
in the condensate. The seeding Bragg pulse created the quasi-particles (in this case condensate excitations in the 
free-particle regime). One can control the momentum of the excited quasi-particles by the angle between the laser 
beams. This could be used to excite phonon-like quasiparticles Jl"5| , and their lifetimes could be determined with the 
pump-probe scheme presented here. 




FIG. 21. Pump-probe spectroscopy of a matter wave grating inside the condensate, (a) First, approximately 5 % of atoms 
were transferred to the recoil state by the two-photon Bragg transition. Then the dynamics of the matter wave grating was 
observed in situ by illuminating the grating with off-resonant dressing light and monitoring the diffracted light intensity. All 
traces were normalized to the same diffracted light intensity at t = 0. The dressing beam intensity was 2.9 mW/cm 2 (bottom), 
5.7 mW/cm 2 (middle), 13mW/cm 2 (top). (b) The initial growth rate of the grating vs. light intensity shows the threshold for 
net gain. The intensity of the dressing beam is given in units of the single-atom Rayleigh scattering rate. Figure is taken from 
Ref. @. 



E. Single-atom and collective behavior 

The optical gain studied above clearly showed the transition from single-atom gain (the first term in the expansion 
in Eq. ^) to collective gain. Varying the intensities of probe and dressing light allows for the study of different physical 
regimes. At low dressing light intensity, below the superradiant threshold, one encounters single-atom behavior, at 
high intensity the system shows collective superradiance. The probe laser intensity determines whether the system 
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shows oscillatory or steady state response, as derived above using optical Bloch equations. Fig. summarizes the 
different regimes. 
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FIG. 22. Schematic diagram of the different regimes of a dressed condensate 
and the probe beams, the dressed condensate occupies different physical regimes 
separated by the threshold to superradiance (G = 1*2, Eq. ^). 
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Probe light traces showing the transition from Rabi oscillations to superradiance are presented in Fig. |23|. As a 
function of the dressing light intensity, the damped Rabi oscillations become faster and almost suddenly turn into a 
giant superradiant pulse. 




i 1 i 1 1 

0.5 1.0 0.5 1.0 

Time [ms] Time [ms] 



FIG. 23. From single-atom Rabi oscillations to collective superradiance. a) The dressing beam intensity was increased from 
2.7, 3.8, 6.1, 9.1, 14.1, 21.3, 32.7, 49.5 to 76 mW/cm 2 (bottom to top). The probe beam intensity was kept at 5 mW/cm 2 . 
b) Numerical solution of the nonlinear optical Bloch equations (Eqs. (78|, ^ and ^). Plotted is for the same experimental 
parameters as in a). 



Previously, recoil related gain based on single-atom phenomena (Recoil Induced Resonances) was observed in 
cold cesium atoms plfl . Collective gain due to the formation of a density grating was discussed as a possible gain 
mechanism for lasing action J5E ] (named CARL — Coherent Atomic Recoil Laser) and pursued experimentally |37 58 
with ambiguous results (see |59 and the discussion in [^0 61 ). Our experiments clearly identify the two regimes and 
their relationship. 

The dressed condensate is a clean, model system for discussing optical and atom-optical properties. The observed 
slow group velocity of the probe laser pulse can be directly related to the dynamics of the amplification process. The 
optical amplification can be described as a reflection of the dressing light by a matter wave grating. The initial delay 
time in the amplification of optical pulses is the time necessary to build up the (quasi-)steady state matter wave 
grating. When the input pulse is switched off, the matter wave grating still exists and diffracts the pump light (as 
observed in Fig. [2l]) creating the trailing edge of the transmitted light pulse. Thus, the slow speed of light is simply 
related to the slow build-up and decay of quasi-particles which we were able to monitor directly. In this microscopic 
picture, all photons propagate with the vacuum speed of light c, the slow group velocity is only a phenomcnological 
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description of the center-of-mass propagation of the amplified pulse. This description leads to the same number of 
photons inside the condensate. Slow light pulses are compressed by a factor c/v g , but the electric field strength is 
unchanged |j^|33). Therefore, the product of the total number of photons within the pulse and the transit time is 
constant. 

Recent demonstrations of slow group velocities for light focused on electromagnetically induced transparency (EIT) 
in a three-level A system [ p0[ . This system features a narrow dip in a broad absorption feature. In our system, the 
broad absorption line is missing. Since the propagation of resonant laser pulses is mainly determined by the narrow 
feature (which determines dni/dw), both systems show analogous behavior (see Fig. |24|). Indeed, if one would add 
a broad-band absorber to the narrow-band amplifier, one would create the same index of refraction as in the EIT 
scheme. 

Although both schemes involve three levels in a A configuration (Fig. ^4|), there are major differences. The am- 
plification scheme does not have a dark state because it has off-resonant couplings to other momentum states which 
are indicated in Fig. ^4|. In the amplification scheme the strong pump pulse connects the initially populated state to 
the excited state in a far-off resonant way. In the EIT scheme the strong coupling laser drives the other leg of the A 
transition and is on resonance. 
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FIG. 24. Comparison between two different methods to generate slow light, (a) Narrow-band optical amplification using a 
two-photon resonance, (b) Electromagnetically induced transparency. 



Finally, we want to resume the discussion of Sect. VI, where we presented the two limits of four- wave mixing of 
light and atoms which correspond to diffraction of light off an atomic grating (for N q > n k - q ) and diffraction of atoms 
off an optical grating (for N q < n k -q). These two limits of atomic stimulation and optical stimulation reflect whether 
the final state occupation number is larger for atoms or photons. Therefore, we have to address what are the relevant 
modes and their occupation numbers. So far, we have disregarded the absolute number of photons. The coupling 
constant C (Eq. p6[ ) is inversely proportional to the mode volume V, hence for the physical description it was only 
relevant to look at the photon density which is proportional to the electric field squared. Therefore, the system can 
be well described semi-classically by the coupling of an electric field (or Rabi frequency) to the atomic density, and 
on the previous pages, we have presented several formalisms for that. 

The question of the photon number becomes relevant when we want to use the concepts of Sect. | where we have 
developed a formalism to describe amplification in a BEC. We showed that the scattering rate W for a condensate 
illuminated by a laser beam is proportional to 

W cx N n k (N q + n k - q + 1). (82) 

This description focuses on two modes for the photons. However, since the experiment takes place in open space with 
freely propagating laser beams, one has to carefully identify the relevant modes for such a description. Here we show 
that all the different cases of light scattering discussed in these notes, Bragg scattering, Rayleigh scattering, and the 
amplification of light and atoms, can be intuitively described by Eq. |82| . 

For Bragg scattering, the optical Bloch equations yield the steady state scattering rate as W = Nqw/2 — = 
Noil'j l /r2- This is in agreement with Eq. [Tci| when we replace the S function by a (normalized) Lorentzian lineshape 
of FWHM of T2, which has a value of 2/7^2 on resonance. Thus the rate for Bragg scattering can be written as 
(neglecting numerical factors and setting S(q) = 1 for large momentum transfer) 

WW « (\C\/h) 2 N Q n k (n k ^ q + l)/r 2 . (83) 
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For the amplification of atoms, the gain equation (Eq. [|5]) implies 

W at0 m_gain » ( | C\ 2 /K)N n k (N q + l)-J^J2,. (84) 

For the amplification of light, the scattering rate at low intensity of the dressing beam is given by Eq. ||— thus we 
are in the regime of scattering atoms off an optical grating. The number of photons n k - q in the input mode is the 
number of photons in a volume V with the cross section A of the condensate and a length Z mo dc ~ C/T2. This choice of 
length is the coherence length of light with a spectral bandwidth T2- This value for Z mo d e can be justified by another 
argument. The coupled equations (Eq. p0| ) reduce the physical description to one mode of the scattered light. The 
angular mode selection occurs by the cigar shape of the condensate. The Bragg scattering process with a linewidth 
T2 couples to only one longitudinal mode if we assume boundary conditions which impose a mode spacing of e.g. 
assuming non-reflecting walls at a distance ps c/T^- To be consistent, the number of photons in the input laser beam 
has to be determined over the same length. 

Now we want to turn to the transition from light amplification stimulated by light to light amplification stimulated 
by atoms. Using Eq. ^ one obtains (again neglecting numerical factors) 

n k - q w nj i Al modc /RX 2 c. (85) 

Expressing the Rayleigh rate R by the atom gain G (Eq. |5^), R ps GA/NqX 2 , we get 

n k - q « N Cl 2 R l mode /Gc w N Q 2 R /GT 2 . (86) 

The number N q of recoiling atoms is approximately the Bragg scattering rate, Q R /T2 times the coherence time I/T2 
times the number of condensate atoms Nq 

N q « N n 2 R /r 2 . (87) 

Comparing Eqs. ^6| and ^7] we see that around the threshold for superradiance (G = T2), the regime N q -C n k ^ q crosses 
over to N q ps n k - q . In other words, the bosonic stimulation factor during the optical amplification, (N q + n k - q + 1) 
(Eq. |32"| ), is approximately n k - q in the linear (single-atom) regime. Near the threshold to superradiance, the N q term 
becomes significant and diverges at threshold. This statement can be generalized to include superradiance. Whenever 
the number of recoiling atoms N q dominates in Eq. [8^ one enters the regime of nonlinear amplification. Without a 
probe beam, i.e. for n k - q — 0, nonlinear scattering and superradiance set in for N q ps 1. 

The photon counting becomes even more confusing when one tries to extend it to the scattered photons. The 
following remarks are rather tentative and reflect hours of controversial and not fully resolved discussions. The 
occupation number of the adiabatically eliminated photon field in Eq. ^ is n k - q — |c*| 2 ps GN q /T c (with N q = \a\ 2 ). 
Since the gain G is the Rayleigh scattering rate into the solid angle fl q of the mode, we can interpret n k - q as the 
photon number emitted into the mode during the damping time 1/T C of the photon field (and stimulated by N q 
atoms in the final state). The physical process is independent of the damping time T c (and therefore of the absolute 
number n k ^ q of photons). In our experiments, the scattered photons hit the walls of the vacuum chamber (or a photo 
detector) within nanoseconds. 

In the case of atom amplification, all scattered photons were emitted spontaneously (in the sense of the superradiant 
cascade described by Dicke |74|| ) and should not be counted in a bosonic stimulation term. Even if these photons were 
emitted into a cavity, they would stimulate the emitting system only after being reflected back. However, it is precisely 
this photon field (Eq. |6^) which leads to the bosonic stimulation by atoms (Eq. |6^) . The Hamiltonian (Eq. ^?]) couples 
atoms only to light. After eliminating the photon field, the (atomic) bosonic stimulation term represents the coupling 
of all atoms to the same mode of the electromagnetic field, and that's how the atoms stimulate each other. Rate 
equations (Eq. [82) only apply after coherences have damped out. The most rapidly damped mode is therefore no 
longer an independent degree of freedom, but slaved to the slower modes. This suggests that the photon number 
n k - q in Eq. |8^ reflects the external laser field, and not the scattered photons — they appear already in the atomic 
stimulation term N q . All these conceptional problems are avoided (but not solved!) when one sticks to the coupled 
equations between the photon and atom fields. 

If we inject probe light near the threshold to superradiance, we have optical amplification by bosonic stimulation of 
atoms! However, the output light intensity is still proportional to the input intensity. In steady state, the population 
N q , the Rayleigh scattering rate (stimulated by these N q atoms) and therefore the number of photons added to the 
probe light are all proportional to the input light intensity. The system is driven by the probe light input. It is 
amplified by a gain factor which includes the dynamics of the coupled atomic and optical fields. 
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Finally, we want to show how to obtain the Rayleigh scattering rate from Eq. |83|. We assume the final photon 
states to be empty (setting nu-q = 0) and multiply with the number of final states. This number Ap is obtained by 
multiplying the density of states Eq. ||| by the linewidth T1T2 and the full solid angle Ap — (dp/dEdfi)(47r)(7ir2). 
This leads to 

WWeigh = (\C\/h) 2 N n k Ap/T 2 a (\C\ 2 /h)N n k ^^4n, (88) 

which is approximately equal to the Rayleigh scattering rate R (Eq. |3^) . The stimulated Bragg scattering rate (Eq. |83| ) 
dominates over the spontaneous Rayleigh rate when the number of photons n k _ q in the Bragg beam is larger than 
the number Ap of accessible modes. Note that this intuitive result is independent of defining modes and boundary 
conditions — the two-photon linewidth T 2 has cancelled out. Subtleties arise only in the case of coherent emission. In 
this case, one has either to describe coherent superposition states of several modes, or identify modes which reflect 
the linewidth of the transition. 



IX. ENHANCEMENT OF SPONTANEOUS EMISSION IN A BOSE-EINSTEIN CONDENSATE 



In the previous sections we have illustrated the rich physics which is described by the rate equations introduced in 
the first section. They suggested that the relevant matrix element is always S(q), and that enhanced scattering relies 
on the population of the final states with several particles or photons. This and the following section broaden this 



picture. This section, which is a slightly modified version of a recently submitted paper |16|, shows that S(q) is the 
matrix element only for scattering particles or light, but not for spontaneous emission. Finally, Sect. [x| discusses that 
collective enhancement is not restricted to bosonic systems and macroscopic population of a quantum state. 

In the previous sections we have shown that light scattering and particle scattering in a condensate can be sup- 
pressed. This suppression is due to the reduction of density fluctuations in a condensate at long wavelengths (longer 
than the healing length). The dissipation-fluctuation theorem implies that the same suppression can be observed 
in the response of a condensate to an external perturbation. Therefore, one might interpret the suppression as the 
"unwillingness" of a condensate against long-wavelength modulations. 

A similar suppression of the long wavelength response (or reduction of the structure factor S(q)) occurs in degenerate 
Fermi gases due to Pauli blocking. Transfer of momentum comparable to or less than the Fermi momentum is 
suppressed due to occupancy of possible final states. Pauli blocking is a single-particle effect and would occur even 
in an ideal non-interacting Fermi gas. In contrast, the suppression of momentum transfer in a BEC is a genuine 
many-body effect due to the atomic interactions and the collective nature of the condensate. It seems that Fermi seas 
and condensates have similar behavior for different physical reasons. 

In addition to particle and light scattering, we want to discuss now the process of spontaneous emission of an 
excited atom initially at rest in the condensate or the Fermi sea. For the Fermi sea, there is suppression: spontaneous 
emission is impossible for recoil momenta less than the Fermi momentum. Surprisingly, the situation is different for 
a condensate where spontaneous emission is enhanced. 

This effect is most easily derived for a homogeneous condensate at zero temperature consisting of N atoms at a 
density p in a volume V. Due to the atomic interactions described by Eq. [H], two atoms in the zero-momentum state 
are coupled to states with momenta +Kk and — Kk. The ground state wavefunction \g) of a BEC has admixtures of 
pair correlations yielding the structure Q 

\g) = \N, 0, 0) - a\N - 2, 1, 1) + a 2 \N - 4, 2, 2) + . . . , (89) 

where a = 1 — 1/u 2 .. Here | No, N k , N-k) denotes a state with No atoms in the zero-momentum state and N± k 
atoms in states with momentum ±hk. In Eq. |8^ a summation over all momenta ?ik is implicitly assumed. The 
average population of momentum states is given by N(k) — ui — 1 = vi, where u k = cosh^fc, v k — sinh^fc and 



tanh20fe = /j,/(hu>^ + p) as already introduced in Sect. 

To study the effect of the presence of a BEC on spontaneous emission, we consider an excited atom at rest added 
to a BEC of N ground-state atoms. This system is described by an initial state |i) = a\ Q \g), where a\ creates an 
electronically excited atom at rest. We use Fermi's golden rule to obtain the rate for spontaneous emission. The only 
difference to the single-atom spontaneous decay rate T comes from the overlap matrix elements to the final momentum 
state |/), {f\a\ L a e fi\i) where a kL is the creation operator for a free ground state atom with momentum of the photon 
fikj_,. Summing over all final states one arrives at 

7bec = T{g\a kL a\ \g) . (90) 
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Thus, the spontaneous emission rate is proportional to the square of the norm of the state vector |e + ) = a* k \g). 
To calculate the norm of |e + ) explicitly, we transform to Bogoliubov operators (see Sect. |l|) and obtain 

= <e + k + } = »L =l + N(k L ) (91) 

where hk s — Mc s is the momentum of an atom moving at the speed of sound. Enhancement of spontaneous emission 
in a BEC is significant if k s becomes comparable to the wavevector kz of the emitted photon since for small momentum 
transfer u\ = k 2 s jk\. 

For comparison, we briefly summarize the suppression of spontaneous emission and light scattering for a fermionic 
system. A Fermi gas at T = with Fermi momentum fikp is characterized by -ZV(k) = 9(kp — k), i.e., all momentum 
states with k < kp — {'otrp) 1 ^ are occupied. If we add an electronically excited atom at rest to the Fermi sea, its 
spontaneous decay rate is suppressed by a factor 

F^ = l-N(k L )=6{k L -k F ). (92) 




FIG. 25. Modification of spontaneous emission (solid line) and light scattering (dashed line) due to quantum degeneracy. In 
a) we have plotted the enhancement factor for spontaneous emission and the suppression factor for light scattering for a weakly 
interacting Bose-Einstein condensate as a function of the light wavevector fct in units of k s , the wavevector of an atom moving 
at the speed of sound. In b) the suppression factors for spontaneous emission and light scattering in a Fermi gas at T = are 
plotted as a function of fez, in units of the Fermi wavevector fc_F- Figure is taken from Ref. 

When off-resonant light with initial wavevector 1cl is scattered from a filled Fermi sphere into an outgoing wave 
with final wavevector Icl + q, the scattering rate is suppressed by |64|] 



S Fc rmi(q) = J dk JV(k)(l - N(k + q)) 

if < q < 2k F , 



(93) 



1 



16fc' 



if q > 2fcp- 



Eq. ^3|is the static structure factor for a Fermi gas at zero temperature. Integrating over all possible scattering angles 
and accounting for the dipolar emission pattern, we find that the total suppression factor for Rayleigh scattering from 
a Fermi sea is given by 
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*a = <! fr2 ' t6 (94) 
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Fig. ^5] shows the influence of quantum degeneracy on the atom-light interaction. Using Eqs. 91, [l^, |9^, and |94| 
we have plotted the rates for spontaneous emission (solid lines) and light scattering (dashed lines), normalized by 
the single-atom rates, for a weakly interacting BEC (Fig. ^5^,) and a degenerate Fermi gas (Fig.p5|b). A significant 
deviation from the free-particle rate is clearly observable if the photon-momentum is comparable to k s for bosons and 
kp for fermions. 

The enhancement of spontaneous emission in a BEC can be intuitively explained as bosonic stimulation by the 
quantum depletion since the enhancement factor (1 + N{k[f)) has the same form as if N{kif) atoms would occupy the 
final state. This intuitive argument is correct, but it would incorrectly predict that light scattering is also enhanced in 
contrast to what we have shown in Sect. H. The suppression of light scattering occurs due to the correlation between 
the admixtures of states with momentum fi,k and — ?ik. This leads to a destructive quantum interference between 
the two processes \N, 0,0) + hq — ► \N, 1,0) and \N — 2,1,1) + Tiq — > | JV, 1,0), in which either an excitation with 
momentum hq is created or an excitation with momentum — hq is annihilated. Both processes transfer momentum 
hq to the condensate and are individually enhanced by bosonic stimulation. Therefore, a simple rate equation model 
would predict enhanced light scattering. However, since the initial states are correlated the two processes leading to 
the same final state interfere destructively for a BEC with repulsive interactions and light scattering is suppressed. 

The static structure factor S(q) of a condensate and the reduced light scattering can be obtained from the Gross- 
Pitaevskii equation as the response to a periodic perturbation pp] |. However, the enhanced spontaneous emission is 
related to the injection of an additional atom into a condensate and appears to be physics beyond the Gross-Pitaevskii 
equation. 

How strong would the enhancement of spontaneous emission in currently realized Bose-Einstein condensates be? 
Condensates of 23 Na atoms confined in an optical trap have reached a density of 3 x 10 15 cm" 3 fl65fl . For this density 
the speed of sound hk s /M — 2.8cm/s and the recoil velocity hk^/M = 2.9cm/s are approximately equal and we 
find N(k]f) 0.15. Thus, the observation of enhanced spontaneous emission in a BEC is within experimental reach. 
Excited atoms at rest could be produced by injecting ground-state atoms with momentum Kki, into a condensate 
and using a counter- propagating laser beam to excite them and bring them to rest. The enhancement of spontaneous 
emission could then be observed as frequency broadening of the absorption line. 

The fact that light scattering is suppressed, but spontaneous emission is enhanced, could be exploited for studies of 
decoherence in a BEC. When a photon is absorbed by a BEC (the first step of light scattering), it creates a (virtual) 
excited state that has an external wavefunction which includes pair correlations. Any decoherence of this coherent 
superposition state, for example by interaction with the thermal cloud, could destroy the interference effect discussed 
above and turn the suppression of light scattering into an enhancement. Another possibility of creating an excited 
state atom in a BEC is using Doppler free two-photon excitation, a scheme already used to probe condensates of 
atomic hydrogen on the Is — > 2s transition [ |66[ . In this case, enhancement of spontaneous emission could be observed 
if the excited state lifetime is longer than the coherence time. 



X. DOES MATTER WAVE AMPLIFICATION WORK FOR FERMIONS? 

Introduction. In the previous sections we have discussed several examples of bosonic stimulation for massive parti- 
cles, bosonically enhanced elastic collisions, superradiance of atoms [O, and matter wave amplification jn||37||. This 
and the four wave mixing of atoms |p8t were described as processes which are bosonically stimulated, i.e., their rates 
are proportional to (N q + 1), where N q is the number of identical bosons in the final state. These experimental 
achievements have raised the question whether these processes are inherently connected to bosonic systems. 

At the Cargese summer school, we presented the view that all these processes do not depend on Bose-Einstein 
statistics and would occur for thermal atoms or even for fermions, although with a much shorter coherence time p[ . 
These suggestions have stirred many controversial discussions at the summer school. This section will reconcile the 
different physical descriptions. The central result is that the stimulated processes mentioned above do not rely on 
quantum statistics, but rather on symmetry and coherence. This section is a slightly extended version of a recently 
submitted paper p9| . 

We also address a widespread misconception about bosonic stimulation which regards stimulated scattering as 
being solely due to quantum-statistical enhancement by the final state, i.e., as if the particles in the final state attract 
other identical particles without any other physical effect. We show that the presence of a macroscopically occupied 
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state increases the density fluctuations of the system, and bosonically enhanced scattering is simply the diffraction 
of particles from these density fluctuations. First, we establish the equivalence of bosonically enhanced scattering, 
diffraction and superradiance which will then be applied to fermionic systems. 

Scattering theory. Sect. | presented basic aspects of the theory of scattering of light or particles from an arbitrary 
system. These results simply followed from lowest order perturbation theory (Fermi's Golden Rule). The double 
differential cross-section for scattering can be decomposed into two factors d ^ ^ = (f^) single S(q : oj). The first one is 
the differential cross-section for the scattering by a single particle (e.g. the Rayleigh cross-section for far-off resonant 
light scattering), the second one is the dynamic structure factor (van Hove or scattering function) S(q,u>) which is 
the Fourier transform of the density-density correlation function: S(q, lo) = (l/27r) J dt e lut (p(q, t)jr(q, 0)) where p(q) 
is the Fourier transform of the particle density operator introduced in Sect. |l|. 

For a non-interacting system of bosons, S(q,u>) can be expressed using the single-particle states |z) with energy Ei 
and occupation numbers Ni as 

S(q,u) = S (q)S(uj)+ 

J2\{j\ e iqr \i)\ 2 Ni(Nj + 1)6 [u - (E, - Ek)/h] . (95) 

The factor (Nj + 1) reflects bosonic stimulation by the occupation of the final state. We have split off the elastic term 
So(q) which describes coherent elastic scattering or diffraction and is simply the square of the Fourier transform of 

the density S (q) = \(pHq))\ 2 = |E W<*| |») | 2 . 

A simple example. It is instructive to apply this formalism to a system of non-interacting bosons which has 
macroscopic occupation in two momentum states with momentum ±fi,k. If the initial state is a Fock state | + 
k) +| — ft) , we find that, apart from forward scattering, the dominant term in S(q,oj) is the bosonically enhanced 
scattering between those two (degenerate) states, S(q,ui) = [N 2 8 q fl + N + (N- + l)8 q -2k + N-(N + + l)<J g ,2fc] 8(v) 
where the Kronecker symbol S q ^ p implies q = p within the wavevector resolution « 1/i of a finite volume with 
length I. Alternatively, we can assume the initial state to be a coherent superposition state \i) with the eigenstate 
\i) = c+| + ft) + c_| — ft) and \c±\ 2 = N±/N and N = N + + N_. Now, the dominant contribution to S(q,ui) comes 
from So(q) = N 2 8 q ,o + N 2 \c+\ 2 \c-\ 2 [8 q ^k + 8 qj -2k] which is equivalent to the Fock state case when the difference 
between N± and N± + 1 can be neglected in the limit of large occupation numbers. 

This equivalence between Fock states and coherent superposition states has been extensively discussed in the context 
of two interfering Bose- Einstein condensates [fzo| — fz2|| and also with regard to optical coherences |fr3| . Those papers 
show that, in many situations, a Fock state is equivalent to an ensemble of coherent states with arbitrary phase. 
Experimental interrogation determines the phase and reduces the ensemble to a single coherent state with a phase 
which will vary from experiment to experiment. For large occupation numbers, one can therefore regard the Fock 
state as an initial state which has not yet "declared its phase" , and, in many cases, for the convenience of calculations, 
replace the Fock state by a coherent superposition state with an arbitrarily chosen phase. 

However, on first sight, the physical interpretation is different. In the Fock state formulation, the enhanced scattering 
results from a macroscopic occupation number in a single quantum state, whereas for the coherent superposition state, 
the scattering is simple diffraction by a sinusoidally modulated density distribution with an amplitude proportional 
to iV|c+c_|. This density modulation acts as a diffraction grating for incident light or particles and has a diffraction 
efficiency proportional to the square of the amplitude. Such a density modulation does not require bosonic atoms. 
It can, for example, be imprinted into thermal or fermionic clouds by subjecting them to a suitable optical standing 
wave. The equivalence of these two descriptions points towards one of the major conclusions of this section, namely 
that macroscopic population of bosonic states is not necessary for enhanced scattering. 

The previous discussion assumed scattering between two degenerate momentum states | ± k) . A simple Gallilean 
transformation generalizes this to two arbitrary momentum states \k±) with energies E±. Now the standing wave 
moves with a velocity h(k + + fe_)/2M where M is the mass of the atoms, and the enhanced scattering appears at 
%uj = ±(E + — _E_) instead of at u> = 0. 

Enhancement of fluctuations. The general results of statistical physics presented above emphasize that enhanced 
scattering must be related to enhanced density fluctuations. Therefore, bosonic enhancement of a scattering rate is 
either due to a density modulation (p(q)} (in the coherent superposition description) or due to density fluctuations 
(in the Fock state description) — the latter can be regarded as a density modulation with an unknown phase. This 
relation allows a more intuitive answer to the question why is there bosonic enhancement when two atoms 1 and 2 
collide in the presence of a condensate with iVo atoms. The standard answer would be that the symmetry of the 
wavefunction enhances the scattering rate into the condensate and into some other state 3 by a factor of (No + !)• An 
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equivalent answer is that the condensate interferes with atom 2 (or 1) and creates a density grating with an amplitude 

1/2 

proportional to N ' which diffracts atom 1 (or 2) into state 3. The grating absorbs this momentum transfer by 
transferring the atom in state 2 (or 1) into the condensate. Therefore, bosonic stimulation can be regarded as 
heterodyne amplification of density fluctuations where the condensate acts as the local oscillator. This alternative 
physical picture emphasizes the role of interference in bosonic stimulation. 

Dicke superradiance. We now want to establish the connection between bosonic enhancement and Dicke super- 
radiance. This will formally introduce the enhancement factor (N + 1) for non-bosonic systems. A system with N 
atoms in two states |±) is conveniently described with the formalism introduced by Dicke to discuss superradiance 
in two- level atoms [Q. It should be emphasized that the only assumption in this treatment is that the N atoms 
couple identically to the probe field (the electromagnetic field or some incident particle beam), i.e., that they have 
the same transition frequency and matrix element without any assumption of quantum statistics. For example, in 
magnetic resonance experiments, the Dicke treatment would apply to different atomic species with the same value of 
the magnetic moment. 

It should be emphasized that Dicke superradiance depends only on the symmetry of the emitting system and does 
not depend on the nature of the emitted particle, whether they are bosons or fermions. For example, the Dicke 
treatment would apply to an ensemble of atoms in an autoionizing state which emit electrons. If the ensemble is 
localized within a de Broglie wavelength of the electron, enhanced superradiant emission would occur. The fcrmionic 
nature of the emitted particles is irrelevant in the so-called microscopic regime |52j where less than one emitted particle 
is in the volume of the sample. This is the situation originally discussed by Dicke. Many experiments on superradiant 
emission of light have been done in regimes with many photons per mode which of course would not be possible with 
fermions. 

Dicke treated the two-level atom as a spin 1/2 system and introduced angular momentum quantum numbers. In this 
subspace, a fully symmetric state of N atoms has spin s = N/ 2 and magnetic quantum number m = (N + —N-)/2. The 
squared matrix element for the transition \s, m± 1) — > \s,m) induced by the ladder operator is (s±m + l)(s=p m). 
Expressing this by initial occupation numbers N±, one obtains N±(N T + 1) |75|-|77|] retrieving the formula of bosonic 
enhancement. The transition rates are largest for the N particle state with s = N/2 which is therefore called the 
state of maximum cooperativity. 

Such a system will couple to the probe field in a superradiant way (i.e. with an up to N times enhanced transition 
rate). In the Bloch vector picture, its dynamics is described as the precession of a macroscopic spin vector with length 
s = N/2. This spin vector decays in a time 1/r where T is the total (homogeneous and inhomogeneous) linewidth of 
the transition |+) — > |— ). Collective superradiant behaviour can only be observed at times shorter than 1/r. 

Matter wave gratings and fermions. Dicke's formalism is usually applied to one-photon transitions between internal 
states, but here we use it to discuss scattering, i.e., a two-photon transition between two momentum states \k±). Let's 
first assume that we have an ideal Bose-Einstein condensate in the k = momentum state. Light scattering between 
momentum states k = and k = q has an infinite coherence time for a non-interacting condensate of infinite size 
(Fig. |26|a) . For a thermal (non-degenerate) cloud of atoms with thermal momentum spread hk t h <C Tiq the transition 
for the transfer of momentum Tiq is Doppler broadened by T = Tik^q/m. For times shorter than 1/r the system will 
behave collectively like the Bose condensed system, i.e., a probe beam would induce transitions between the k = 
and k = q momentum states at a rate proportional to Nk—a{Nk— q + 1) where iVfc=o,g refers to the total number of 
atoms in states with momentum around k = 0, q. 

Once we have distributed the particles over many initial states, indistinguishability and quantum statistics don't 
play any role. Therefore, the only modification for a Fermi degenerate cloud is to replace fcth with the Fermi wavevector 
&f in the expression for the inhomogeneous broadening (Fig. p6|b). Due to the assumption Tikp <C Tiq, Pauli blocking 
due to scattering into already occupied states is absent. If this assumption is not made, a part of the cloud becomes 
inactive, and our discussion would apply only to the atoms near the Fermi surface. 

The previous paragraph generalized the bosonic Fock state ensemble to non-degenerate and fermionic clouds. We 
now come back to the coherent superposition state. For bosons, it can be produced from a Bose-Einstein condensate 
in the k = state by applying a (so-called Bragg) pulse of two laser beams which differ in wavevector by q and in 
frequency by the recoil frequency Tiq 2 /2m. Those beams resonantly drive the transition between momentum states 
k = and k = q j78pC| ] and prepare the superposition state discussed above. Similarly, in a thermal (or fermionic) 
cloud, the Bragg pulse creates a modulated density distribution with wavelength 27r/g which has the same contrast 
as in the bosonic case and will diffract light or atoms at the same rate. However, due to the thermal motion with 
velocity hk^/M, this grating decays during a time M/Tik^q = 1/r (for the fermionic case, kp has to be substituted 
for fcth)- Thus the Dicke picture and the diffraction picture agree. 
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FIG. 26. Momentum transfer q (a) to a Bose-Einstein condensate, (b) to a Fermi sea, and (c) to a momentum squeezed 
degenerate Fermi cloud. Shown are the populated states vs. the k- vector. The momentum spread fc_F of the Fermi sea introduces 
Doppler broadening of the transition and a finite coherence time, whereas the coherence time in (a) and (c) is infinite. Figure 
is taken from Ref. [p9|. 



Coherence time. The Doppler broadening discussed above seems to imply a fundamental limit to the coherence 
time of a Fermi system. However, at least in principle, one can prepare a Fermi system with infinite coherence time by 
starting out with a cloud which is in a single momentum state along the z axis, but occupies many momentum states 
along x and y. With a Bragg pulse transferring momentum hqz, one can prepare a system which shows collective 
behavior for scattering particles or light with momentum transfer hqz with an infinite coherence time (Fig. |26|c). 
Therefore, there is no direct connection between a long coherence time and a high phase-space density. In this 
ensemble, the scattering is between the states \k z = 0) ® \k x , k y ) and \k z = q) ® \k x , k y ). Therefore, we have enhanced 
scattering into the \k z = q) quantum state, but the atoms may differ in other quantum numbers. What matters is 
only the symmetrization of the many-body wavefunction along z. The other quantum numbers ensure that there is no 
conflict with the Pauli blocking for fermionic systems. This is analogous to the separation of electronic wavefunctions 
into a symmetric part (e.g. the spin part) and an antisymmetric part (e.g. the spatial part) where the coupling to an 
external field (e.g. electron spin resonance experiment) only depe nds on the symmetric part. 

Experiments. The experiments both on superradiance (Sect. VII and 14|) and four-wave mixing pfifl in Bose-Einstein 
condensates have in common that a matter wave grating formed by two macroscopically occupied momentum states is 
probed, either by light or by atoms. Both experiments create the coherent superposition state discussed above using a 
Bragg pulse. In the limit of low intensity of the probe beam, the scattering is independent of the nature of the probe 
particles — one could have used any kind of radiation, bosons or fermions |7S(] . The bosonic stimulation observed in 
both experiments demonstrates the dynamic nature of the matter wave grating. Each time, a particle or photon is 
diffracted, the amplitude of the grating grows. 

In practice, it is difficult or impossible to carry out these experiments with fermions or thermal atoms. When we 
observed superradiance of a condensate, we couldn't observe similar behaviour above the BEC transition temperature 
since the threshold laser intensity for superradiant gain is several orders of magnitude higher (see Ref. [H for details) . 
Furthermore, the superradiance may be suppressed by heating or other decoherence processes. The shorter coherence 
time for non-BEC samples should be even more crucial for the four-wave mixing experiment where the matter wave 
grating is probed by very slow atoms which have a long transit time of about 1 ms through the sample. Another 
concern are incoherent processes which accompany the stimulated processes discussed so far. Since the incoherent 
processes scale linearly with the number of atoms, whereas the stimulated process is proportional to TV 2 , there is in 
principle always a regime where the stimulated process dominates 

Discussion. Coming back to the initial question: Is matter wave amplification possible for fermions? The answer 
is yes, if the system is prepared in a cooperative state and the amplification is faster than the coherence time. 
However, this amplification does not pile up atoms in a single quantum state, but rather in states which are in the 
same (or approximately the same) momentum state along z, but differ in other quantum numbers. Therefore, this 
amplification can be regarded as amplification of a density modulation or as amplification of spatial bunching of 
atoms. Alternatively, one can regard the density modulation as a collective excitation of the system which involves 
bosonic quasi-particles (e.g. phonons). Superradiance and four- wave mixing (both with bosons and fermions) can 
then be ascribed to bosonic stimulation by those quasi-particles. 



3 In the case of a fermionic sample, the density of particles can only be increased by increasing the momentum spread of the 
sample. This can increase the incoherent "background" for four-wave mixing. For the ensemble in Fig. ^c an increase in N 
at constant volume requires the transverse velocity spread to increase as N 1 ^ 2 . Therefore, for large N, the incoherent elastic 
scattering rate increases as iV 3//2 , still more slowly than the stimulated scattering. 
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The phase-coherent matter wave amplification for fermions would start with a short Bragg pulse which puts some of 
the atoms into a recoil state which is then amplified. This superposition of two momentum states creates a matter wave 
grating. This can be regarded as the interference pattern of each atom with itself with all the individual interference 
patterns being exactly in phase. Matter wave amplification occurs when a single laser beam is diffracted off this 
grating increasing the amplitude of each atom to be in the recoiling state. Therefore, the matter wave amplification 
scheme of Refs. |l7],[37| would work for fermions, provided the whole process can be done in the short coherence time 
of the fermionic matter wave grating. 

Of course, there is a fundamental difference between bosons and fermions which is reflected in the symmetry of 
the total wavefunction. A bosonic system with two macroscopically occupied quantum states is always in a fully 
symmetric and maximally cooperative state. In other words, if two independent Bose condensates cross each other, 
there is always a macroscopic interference pattern (as observed experimentally p3|), which is reflected in S(q, u>) being 
proportional to N 2 (or to N+N-, to be more precise). It is this density modulation which can be amplified by the 
dynamic diffraction discussed in this paper. If two beams of fermions overlap, there is no macroscopic interference, 
unless the two beams were prepared in a symmetric way, e.g. by generating one of the beams by a Bragg pulse from 
the other one. 

Our discussion of scattering without change of the internal state can be generalized. For example, if atoms scatter 
into the condensate through a spinflip process, the density grating has to be replaced by a polarization or coherence 
grating. Such gratings were experimentally studied for laser-cooled atoms 



This discussion has focused on bosonically enhanced scattering. Similarly, bosonic enhancement of spontaneous 
emission depends only on a cooperative initial state and not directly on quantum statistics. For scattering, the 
relevant coupling strength are the density fluctuations. For spontaneous emission, it is the electric dipolc moment. 
Both are enhanced by the presence of a Bose condensate, in the latter case because the excited atom corresponds to 
a Dicke vector of spin s = N/2, m = —(N/2) + (1/2) which couples more strongly to the vacuum fluctuations of the 
electromagnetic field than an individual atom. Alternatively, the enhanced spontaneous emission can be regarded as 
the constructive interference of an "emitted" ground state atom with the macroscopic ground state matter wave. This 
picture is analogous to the semi-classical interpretation of stimulated emission of light. Ref. |75| shows that bosonic 
stimulation of photons is due to the constructive interference of the emission of a classical oscillating dipole with the 
incident field in the forward direction. 



XI. DISCUSSION 



This paper has summarized our recent experiments on Bose-Einstein condensation with the unifying theme of 
enhancement and suppression. Suppression of scattering or dissipation can arise for two different reasons. The phonon 
and vortex nature of the collective excitations together with energy and momentum conservation allow dissipation 
only above a critical velocity. In addition, one has to consider the dynamics of the excitation process. For microscopic 
particles, this is reflected in the matrix element S(q) which characterizes how easily can the condensate absorb 
momentum in a scattering process. For macroscopic motion, it is reflected in a critical velocity for vortex nucleation. 
Scattering processes are also enhanced by the population of the final states (bosonic stimulation) . Optical stimulation 
by a laser beam was used in Bragg scattering, and matter wave stimulation led to superradiance and matter wave 
amplification. 

We have also discussed some subtleties which go beyond the simple picture using rate equations and occupation 
numbers. A condensate in its ground state is in a coherent superposition state of the zero- momentum state with 
correlated pairs with momenta ±q (the quantum depletion). We have shown that the population in the quantum 
depletion can cause bosonic stimulation of spontaneous emission. However, for a scattering situation, there are two 
bosonically enhanced pathways which destructively interfere (causing S(q) < 1). Therefore, the concept of bosonic 
stimulation can be applied to the many-body state of a condensate, but with caution. Finally, we have shown that 
some form of matter wave amplification is possible in fermionic samples. This required a careful discussion of bosonic 
stimulation by particles vs. quasi-particles and the role of symmetry vs. quantum degeneracy. 

In closing we want to point out that the rich and complex physics displayed here is based essentially on two four-wave 
mixing Hamiltonians. One describes s-wave interaction of the condensate which is responsible for all the many-body 
effects discussed here including superfluidity. The other one couples the atoms to the light (or impurity atoms) and 
led to superradiance and amplification of light and atoms. 
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